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ABSTRACT 


This investigation studies the short-wavelength buckling (or 
the microbuckling) and the interlaminar and inplane shear failures of 
multi-directional composite laminates loaded in uniaxial compression. 

A laminate model is presented that idealizes each lamina. The fibers 
in the lamina are modeled as a plate, and the matrix in the lamina is 
modeled as an elastic foundation. The out-of-plane w displacement 
for each plate is expressed as a trigonometric series in the half- 
wavelength of the mode shape for laminate short-wavelength buckling. 
Nonlinear strain-displacement relations are used. The model is 
applied to symmetric laminates having linear material behavior. The 
laminates are loaded in uniform end shortening and are simply 
supported. 

A linear analysis is used to determine the laminate stress, 
strain, and mode shape when short-wavelength buckling occurs. The 
equations for the laminate compressive stress at short-wavelength 
buckling are dominated by matrix contributions. The effects of fiber 
volume fraction on the compressive stress at short-wavelength buckling 
is reported for a laminate with any stacking sequence and any 


xv 



xvi 


thickness. The laminate mode shape at short-wavelength buckling is 
discussed. 

A nonlinear analysis for laminae with initial imperfections is 

used to determine laminate stresses and interlaminar strains. Results 

are presented for imperfection-amplitude-to-lamina-thickness ratios of 

0.1 and 0.5. The nonlinear behavior of [ 0_ ] , [0/90] , [±-45 ] , and 

s s s 

[ +45/0/— -45/90 ] s laminates is discussed. Results are presented for 

laminate stress, end shortening, and w displacement. The w 

displacement gradients cause significant interlaminar shear strains 

Y at laminate stresses that are much lower than the laminate stress 
xz 

at short-wavelength buckling, and the distribution of these strains is 

described. The effect of fiber volume fraction on the nonlinear 

laminate response is presented for [0_] and [±^5] laminates. 

c. s s 

A failure criterion for compression-loaded laminates is 

discussed. Laminate failure that is initiated by outer-lamina 

buckling, by interlaminar shear strains from lamina imperfections, or 

by inplane matrix shearing is considered. The laminate strength is 

calculated as a function of lamina orientation for [ +0 ] laminates. A 

s 

simple method referred to as the stiffness-ratio method is described 

O 

for predicting the strength of 0 -dominated laminates. 



CHAPTER 1 


INTRODUCTION 


Composite materials are being widely used in a variety of 
applications ranging from sporting equipment to primary structures for 
commercial transport aircraft. The design of composite structures is 
often a combination of traditional practices used for metal structures 
and empirical criteria. This design technique is easy to use but 
neglects the unique mechanisms that dominate the behavior of composite 
structures. The most efficient composite structures are designed 
using a thorough understanding of the basic response mechanisms of 
these materials. 

The uniaxial compressive strength of a composite laminate is a 
fundamental property. The mechanics of uniaxial compressive failure 
have been studied by many researchers for unidirectional laminates 
having fibers parallel to the loading direction. These laminates fail 
when short-wavelength buckling (or microbuckling) of the fibers occurs 
or when shear failure at the fiber-matrix interface occurs. Few 
similar investigations have been conducted for multi-directional 
composite laminates although multi-directional laminates are more 
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widely used than unidirectional laminates in composite structural 
components. This study investigates the short-wavelength buckling of 
multi-directional composite laminates loaded in compression. This 
study also addresses the interlaminar shear failures due to short- 
wavelength imperfections and the inplane shear failures in these 
laminates. The results of this study are reported herein. 

1 .1 Chronological Review of the Literature 


This section reviews the literature on models for predicting 
the compressive strength of composite materials. This review focuses 
on studies that emphasize the short-wavelength buckling phenomenon and 
the interlaminar shear failures due to short-wavelength initial 
imperfections. The review is chronological except that similar 
subsequent publications by an author are discussed with the first 
reference. A summary of the review is found at the end of this 
section. 

In I960 Dow and Gruntfest [1] postulated that the compressive 
failure of unidirectional laminates was the result of either of the 
following two phenomena: (1) high tensile stresses perpendicular to 

the loading direction and at the fiber-matrix interface; or (2) 
buckling of the fibers within the matrix. The former phenomenon is 
often called transverse tension failure, and the latter phenomenon is 
often called microbuckling. An equation to predict the laminate 
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compressive strength at microbuckling was derived using reference 2 
and is based on a model consisting of columns supported by an elastic 
foundation. An empirical constant was included in this equation. 

This reference was the first to associate fiber instability with the 

i 

compressive strength of unidirectional laminates. 

Fried [ 3 ], Fried and Kaminetsky [4], and Fried [5] studied the 
influence of the matrix, the reinforcement, and the matrix- 
reinforcement interface on the compressive failure of composite 
materials. Experimental results from unidirectional laminates were 
reported in reference 3 and led Fried to suggest the following failure 
sequence: the reinforcement phase in a composite material carries the 

compressive load until the rigid matrix phase yields; upon yielding, 
the matrix flows and no longer supports the reinforcement; the 
reinforcement buckles; and the composite material fails 
catastrophically. The laminate compressive strength was expressed as 
a linear function of the matrix yield stress by assuming that (1) the 
maximum stress in the matrix is its yield stress and that (2) the 
strain in the matrix is the same as the strain in the reinforcement. 

In his investigation Fried also identified the straightness of the 
reinforcement and the bond between the reinforcement and the matrix as 
secondary effects on the compressive strength. 

Experimental results in a subsequent study [A] supported 
Fried's initial conclusions. The results in reference 4 showed that 
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the compressive strength for unidirectional' laminates increases with 
increasing matrix yield stress for the materials studied. The yield 
stress for each matrix material was presented, but no other properties 
(e.g., Young's modulus or shear modulus) of the matrix were 
discussed. 

Fried [5] also studied the compressive failure of filament 

wound [0/90] -class laminates. He observed interlaminar cracking in 
s ° 

failed compression specimens. This failure mode for specimens loaded 

in compression was similar to the failure mode for specimens loaded in 

interlaminar shear, and Fried assumed that failure was the result of 

debonding at the matrix-reinforcement interface for both load cases. 

Fried suggested that debonding initiated at voids. Experimental 

results showed that laminate compressive strength was inversely 

proportional to void content. As discussed in references 3 and 4, the 

compressive strength for laminates with very low void contents was 

limited by the matrix yield stress. 

Leventz [6] studied the influence of fiber diameter on the 
compressive strength of glass-epoxy composite materials. Experimental 
results showed that the compressive strength was maximized using a 
0.005-in. -diameter glass fiber. The range of fiber diameters for this 
study was from 0.001 inches to 0.010 inches. Leventz reported that 
the failure mode for the unidirectional composite materials in his 
study appeared to be a combination of fiber instability and matrix 
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shear failure. He suggested that the laminate compressive strength 
was proportional to the square of the fiber diameter although the 
experimental results did not support this theory. Leventz reported 
that the larger diameter fibers were more collimated than the smaller 
diameter fibers, but he did not investigate the influence of fiber 
straightness on the laminate compressive strength. 

In 1965 Rosen [7] reported on his classic study of the 
compressive failure of unidirectional composite materials. This study 
focused on fiber instability in glass-epoxy laminates using a two- 
dimensional model. The fibers were modeled as columns supported by an 
elastic matrix foundation. Rosen suggested that compressive failure 
for a unidirectional composite material occurred when the fibers 
buckled into either of two possible short-wavelength modes, the 
extension mode or the shear mode. These mode shapes are shown in 
figure 1.1. For the extension mode, the deformation of the matrix 
material between fibers is extension in the direction perpendicular to 
the fibers. For the shear mode, shear deformation occurs in the 
matrix material. The wavelength is short compared to the length or 
the width of the test specimen for both buckling modes. Using an 
energy formulation, Rosen obtained 


2V, 


V. E m E r 
f m f 

3(1 -V f ) 


1/2 


( 1 . 1 ) 
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for the extension mode and 


m 


1-V, 


E f h 


12A 


( 1 . 2 ) 


for the shear mode where 

o Q = compressive strength of the unidirectional composite 
material 

V f = fiber volume fraction 

E = Young's modulus for the matrix 
m 

E^. = Young's modulus for the fiber 

G = shear modulus for the matrix 
m 

h = fiber diameter 

A = half-wavelength of the buckling mode shape 

The magnitude of the second term in equation (1.2) is small compare', 
to the magnitude of the first term since the half-wavelength is much 
larger than the fiber diameter. Neglecting the second term, Rosen 
obtained the approximate (and more familiar) equation for the shear 
mode 


a 


c 


G 


(1.3) 
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The composite material compressive strength was predicted using the 
lower value from equations (1.1) and (1.3). Equation (1.3) gives the 
lower prediction for most composite materials. 

Rosen recognized that equation (1.3) gave predicted strengths 
that were two to three times greater than experimental strengths when 
the elastic shear modulus of the matrix was used. He obtained more 
realistic strength predictions by assuming the matrix shear modulus 
was a function of the applied load. 

Unlike previous theories, Rosen's analysis did not use 
empirical factors for predicting compressive strength. Also, equation 
(1.3) illustrates that the composite material compressive strength is 
a function of the matrix shear modulus for fiber-instability-initiated 
failures. 

The results from reference 7 were also included in a 
subsequent report by Dow and Rosen [8]. Also, Dow, Rosen, and Hashin 
[9] modified equation (1.3) to account for an elastic-perfectly- 
plastic response of the matrix. The matrix material was assumed to be 
isotropic and incompressible after yielding. The authors of reference 
9 obtained 


a c = 


V f E f o f 


1/2 


3(1-V f ) 


(1.4) 
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where is the yield stress of the matrix, and the other variables 

are defined above. The use of equation (1.4) is limited since many 
matrix materials do not exhibit elastic-perfectly-plastie behavior. 

Ekvall [10] also studied fiber instability for glass-epoxy 
laminates loaded in compression. He observed that the matrix yielded 
well before laminate failure. This observation contradicted the 
theory in reference 3. Ekvall predicted the buckling wavelength for 
fiber instability to be on the order of a fiber diameter, and this 
prediction agreed with his experimental results. 

Hayashi [11] reported on a "shear instability failure" for 
orthotropic materials. He postulated that such a failure occurred 
when a material's flexural rigidity was significantly greater than its 
shear rigidity. This concept of a shear-deformation-dominated 
instability for orthotropic materials is the same as that discussed by 
Rosen [7] for fiber-reinforced composite materials. Hayashi minimized 
the potential energy of a compression-loaded orthotropic material and 
obtained an expression for the strength associated with the shear 
instability failure 


= G 


(1 .5) 
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where o Q is the compressive strength and' G is the shear modulus of 

the material. The shear modulus for a unidirectional composite 
material can be calculated as a function of constituent properties 
from 


G 

c 



G 


m 

+ d-v f ) 


(1 . 6 ) 


where G c is the shear modulus for a unidirectional composite 

material and G^, is the shear modulus of the fiber (see reference 

12). The other variables are defined in the preceeding paragraphs. 

The shear modulus of the fiber is much greater than the shear modulus 
of the matrix for most fiber-reinforced composite materials, and the 
first term in the denominator of equation (1.6) can be neglected. The 
shear modulus for a unidirectional composite material is approximated 
by 


G =* 
c 



(1.7) 


When equation (1.7) is substituted into equation (1.5) (Note: G = G c 
for composite materials), equation (1.3) is obtained. 
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Independent of Rosen, Schuereh [13] also obtained equations 
(1.1) and (1.3) for the short-wavelength buckling of unidirectional 
composite materials. Schuereh tested boron-magnesium laminates and 
obtained good correlation between theory and experiment. 

In 1966 Foye [14] showed that the upper bound for any 
material's compressive strength was its shear modulus (ef., equation 
(1.5)). Also, he showed that this upper bound could be approximated 
by equation (1.3) for unidirectional composite materials. Foye 
suggested that differences between experimental strengths and 
predicted strengths for composite materials were due to "local 
imperfections." He observed fiber instability in randomly located 
regions throughout the material. He also modified his compressive 
strength predictions by including the effects of voids and matrix 
fillers. 


Some researchers attempted to isolate fiber instability in a 
laminate by studying a single fiber surrounded by matrix. Hermann, 
Mason, and Chan [15] developed a beam model that included initial 
waviness of the beam. The authors of reference 15 found that the 
short-wavelength buckling loads predicted by their model agreed with 
similar loads for the extension mode in reference 13. Also, the 
response of the model was significantly influenced by the initial 
waviness of the beam. Sadowsky, Pu, and Hussain [16] developed a 
model to study fiber instability caused by manufacturing-induced 
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residual thermal stresses. The authors of reference 16 predicted the 
compressive strain for short-wavelength buckling of the fiber. They 
showed that the shear deformation of the fiber had a negligible effect 
on the fiber compressive strain at buckling for most composite 
materials. The predicted strains appear to be more than an order of 
magnitude greater than experimental results. 

Crawford [17] studied the compressive behavior of a boron- 
polyimide-f ilm layered composite material. Analytical and 
experimental results were presented. Crawford modified the analysis 
of reference 13 to account for any inplane-load-carrying capability of 
the matrix. The inplane axial stiffness of the boron layer was 
expressed as a function of the layer's initial waviness. Crawford 
suggested that this initial waviness caused interlaminar normal and 
shear stresses that resulted in laminate failure prior to short- 
wavelength buckling. He predicted the amplitude and half-wavelength 
of the initial waviness for this laminate to be approximately five 
times and eighty times the thickness of the boron layer, 
respectively. 

Yue, et al. [18] studied the compressive behavior of A1 -CuA1 2 

eutectic composite materials. The authors of reference 18 expanded 
the results of reference 9 to model linear strain-hardened matrix 
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materials and obtained reasonable agreement between theoretical and 
experimental strengths. 

Lager and June [19] used Rosen's results [7] to experimentally 
study the effect of matrix modulus and of fiber volume fraction on the 
unidirectional compressive strength of boron-epoxy laminates. Two 
matrix materials were used to fabricate specimens with fiber volume 
fractions ranging from 0.05 to 0.46. The authors of reference 19 
modified equations (1.1) and (1.3) by including an "influence 
coefficient" that was determined empirically. They used the influence 
coefficient to account for softening of the matrix with increasing 
applied load. The coefficient was the same for both equations. The 
strength of specimens fabricated using the stiffer matrix were higher 
than the strength of similar specimens fabricated using the softer 
matrix. The experimental strengths for specimens having fiber volume 
fractions less than ten percent agreed with the predicted strengths 
from modified equation (1.1), and the experimental strengths for 
specimens having fiber volume fractions greater than ten percent 
agreed with the predicted strengths from modified equation (1.3). 

Chung and Testa [20] applied Biot's mechanics of incremental 
deformations [21] to a model for fibrous composite materials in their 
study of short-wavelength buckling. The model consisted of fiber 
beams supported by an elastic matrix foundation. The authors of 


reference 20 calculated the laminate stress that buckled fibers into 



either the extension mode or the shear mode (figure 1.1). The 
equations for these stresses reduced to Rosen’s equations [7] (i.e., 
equations (1.1) and (1.3)) when the buckling wavelength was much 
greater than the thickness of the matrix between fibers. The authors 
of reference 20 also presented experimental results from specimens 
with fiber volume fractions less than 0.10 that agreed with their 
analytical results. The equations in reference 20 are useful but much 
more difficult to use than the equations in reference 7. 

Russian researchers have also predicted the compressive 
strength of unidirectional composite materials [22-27]. They obtained 
equation (1.5) using nonlinear elasticity theory [22-24,26,27]. 

Skudra, Kalnays, and Bulavs [25] postulated that the laminate 
compressive strength was strictly a function of compressive and shear 
failures of the composite material's constituents. Limited 
experimental results are also presented in references 24 and 25. 

De Ferran and Harris [28] studied the compressive strength of 
steel-wire - polyester-resin composite materials. They suggested that 
the laminate compressive failure was a function of the wire (cr fiber) 
tensile strength and that a simple ”rule-of-mixtures" analysis could 
predict the compressive strength. They obtained experimental 
strengths that agreed with their predicted strengths although the 
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analysis neglects any of the fiber instabilities or matrix failures 
documented in previous studies. 

Kiusalaas and Jaunzemis [29] also used incremental deformation 
theory [21] to study short-wavelength buckling of composite materials. 
Their investigation was motivated by the beam-on-an-elastic-foundation 
model used by Chung and Testa [20]. Kiusalaas and Jaunzemis derived a 
continuum theory for buckling of a laminated medium. The medium 
consisted of alternating fiber and matrix layers. The results of 
reference 29 reduce to the results of reference 20. 

Hayashi [30] and Hayashi and Koyama [31] expanded Hayashi's 
previous theory for the compressive strength of orthotropic materials 
[11] to unidirectional composite materials. Hayashi [30] assumed that 
a compression-loaded composite material fails when the axial stress ir 

the matrix equals the "shear instability limit" of the matrix, o . 

He defined this limit as the shear modulus of the matrix and 
calculated the compressive strength of a unidirectional composite 
material by assuming equal axial strain in the fiber and the matrix. 
Hayashi obtained the rule-of-mixtures equation 


o = a r V + a ( 1-V.) 
c f f m I 


( 1 . 8 ) 
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where is the fiber stress corresponding to the axial strain at 

failure. The other variables are defined in the preceeding 
paragraphs. The fiber had linear elastic material response. This 
analysis considers matrix instability (instead of fiber instability) 
as the cause of laminate failure. The predicted strengths from 
equation (1.8) were compared to experimental strengths in references 
30 and 31 . Only some of the data agreed with the predictions, and no 
explanation was given for discrepancies between analytical and 
experimental results. 

Independent of Hayashi, Orringer [32] also postulated that a 
matrix instability could cause compressive failure of composite 
materials. Orringer defined a "matrix critical strain" and calculated 
the composite material's strength at this strain using a rule-of- 
mixtures equation. He also studied short-wavelength buckling of 
fibers. Orringer suggested that "progressive buckling" of the fibers 
led to laminate failure. He assumed that the progressive buckling 
started at voids or at disbonds within the laminate. 

In 1972 Greszczuk [33, 3^1 reported on his experimental 
investigations of the effects of the constituents on the compressive 
strength of fiber- and lamina-reinforced composite materials. These 
investigations include the most thorough experimental results to date. 
The author varied the matrix material, the fiber material, the fiber 
diameter, the fiber volume fraction of the specimen, the initial 
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imperfections (i.e., initial waviness) in the specimen, the specimen 
geometry, and the specimen's boundary conditions for unidirectional 
fiber-reinforced specimens. Greszczuk performed similar studies for 
lamina-reinforced specimens with 6061 -T6 aluminum laminae. He 
reported on the influence of constituent properties on the composite 
material's compressive failure mode. Short-wavelength buckling of the 
reinforcement occurred when the shear modulus of the matrix was low. 
Compressive strength failure of the reinforcement occurred when the 
shear modulus of the matrix was high. Initial imperfections 
significantly decreased the strength of all specimens. The results 
presented in references 33 and 34 were also presented in subsequent 
reports [35-391. 

Suarez, Whiteside, and Hadcock [40] studied the compressive 
strength of multi-directional boron-epoxy laminates. This study was 
one of the first to investigate lamina instability in multi- 
directional laminates. The authors of reference 40 assumed that the 
outer-most laminae buckled in a short-wavelength mode; the remainder 
of the laminate was treated as an elastic foundation with both 
extensional and shearing stiffness. The fibers in the outer-most 
laminae were always oriented parallel to the direction of loading 
(i.e., 0-degree laminae). Some experimental results agreed with 
predicted strengths when the initial waviness of the outer laminae was 
included in the analysis. The authors assumed that the ratio of 
waviness amplitude to lamina thickness ranged from 0.2 to 0.4. 



17 


Filament failure and global instability of the specimen were also 
discussed. 

Hackett, Tarpy, and Wood [41] studied fiber instability by 
performing tests on a single steel wire embedded in a block of epoxy 
resin. They obtained the load and stress distribution at buckling 
using a photoelastic stress analysis and compared these experimental 
results with analytical results they obtained using a finite element 
analysis. The authors of reference 41 reported reasonable agreement 
between the experimental and the analytical results. This study 
illustrates a useful test technique but has limited application to the 
short-wavelength buckling of composite laminates. 

Lanir and Fung [42] considered the buckling and postbuckling 
response of cylindrical columns of matrix reinforced with parallel, 
straight fibers. They suggested that the compressive failure of a 
unidirectional laminate occurred at a load much greater than the 
buckling load of the fiber. They did not consider a short-wavelength 
buckling mode as the initial mode shape for the buckled fiber. The 
buckling mode shape changed with increasing load (the wavelength 
decreased with increasing load). They assumed that failure resulted 
from high stresses at the fiber-matrix interface. 

Independent of references 11 and 14, Kao and Pipkin [43] 
showed that the critical stress for short-wavelength buckling was 
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equal to the shear modulus of the composite material (see equation 
(1.5)). They also showed that the shear mode for short-wavelength 
buckling was an admissible deformation for fiber-reinforced columns. 
This study discusses short-wavelength buckling as a basic response 
phenomenon for unidirectional laminates; the compressive failure of 
such laminates is not discussed. 

Kulkarni, Rice, and Rosen [44,45] modified equation (1.3) to 
account for shear deformation of the fiber and for imperfect bonding 
at the fiber-matrix interface. They obtained an expression for the 
compressive strength of a unidirectional laminate, 




1 -d-k)v 




i 




(1 .9) 


where is the shear modulus of the fiber and k is a bonding 

parameter. They defined the range for k as 

| $ k < 1 (1.10) 



where the limits represent no bonding and perfect bonding, 
respectively. The authors showed that the compressive strength 
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predicted by equation (1.9) was less than the compressive strength 
predicted by equation (1.3). 

Davis [46,47] included both nonlinear material behavior and 
initial waviness of the laminae in his analysis of the compressive 
strength of unidirectional composite materials. A model consisting of 
fiber layers and matrix layers was analyzed as a multilayered 
Timoshenko beam loaded in compression. The matrix material had a 
nonlinear shear stress-strain behavior. Davis measured the initial 
waviness of the laminae in his boron-epoxy specimens, and the ratio of 
initial waviness amplitude to lamina thickness ranged from 0.2 to 1.6. 
He assumed that failure was caused either by delamination or by short- 
wavelength of the lamina. Davis referred to thes short-wavelength 
buckling as a shear instability. He found that shear instability was 
the dominant cause of failure for these specimens. Compressive 
strength predictions based on equation (1.5) agreed with the 
experimental results. 

Hanasaki and Hasegawa [48] and Wang [49] independently 
conducted studies similar to that conducted by Davis [46]. The 
results of references 48 and 49 are also similar to the results of 
reference 46. 

Evans and Adler [50] suggested that fiber kinking was the 
dominant failure mechanism for unidirectional laminates. Kinking 
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occurred subsequent to short-wavelength buckling within the laminate. 
The authors of reference 50 studied the mechanics of fiber kinking and 
obtained a approximate expression for the "critical kink formation 
stress." They presented experimental results from three-dimensional 
carbon-carbon woven composite material specimens to illustrate kink 
geometry. 


Similar to Evans and Adler, Maewal [51] studied the short- 
wavelength buckling and postbuckling behavior of unidirectional 
laminates in compression. He used a three-layer model consisting of a 
fiber layer surrounded by matrix layers to analyze the shear mode 
buckling and initial postbuckling of these laminates. Maewal found 
that his model predicted a postbuckling stiffness for the laminate 
that was approximately one-third of the prebuckling stiffness of the 
laminate. He suggested that initial waviness of the fiber layer did 
not significantly affect the short-wavelength buckling stress. These 
results contradict previous studies that show unidirectional laminates 
have no postbuckling stiffness and show that initial waviness of the 
laminae significantly affects the stress distribution in the laminate 
(e.g., reference 46). 

Budiansky [52] used the results of the previous two studies 
[50,51] in his investigation of fiber kinking. He included matrix 
plasticity and initial waviness of the fibers in his analysis of 
kinking. Budiansky identified the shear stiffness and shear strength 



21 


of the composite material as the most important parameters affecting 
"kink strength." This reference contains a good discussion of the 
importance of understanding the role of the matrix and the role of 
initial waviness of the fibers in predicting the compressive strength 
of unidirectional laminates. 

Hahn and Williams [53] studied the compressive failure 
mechanisms of unidirectional composite materials using a nonlinear 
model that included initial waviness of the fiber (called fiber 
curvature) and nonlinear material behavior for the matrix. The 
authors of reference 53 determined the loading on a single fiber 
surrounded by matrix and obtained 


for the compressive strength of a unidirectional laminate with a 
linear elastic matrix and straight fibers. In equation (1.11) G c is 

the composite material shear modulus. A laminate fails when short- 
wavelength buckling occurs. Equation (1.11) differs from previous 
predictions [7,11] by the factor V f . Hahn and Williams suggest that 

this difference is the result of the free body diagram used in the 
analysis. As expected, they found that the compressive strength of 
composite laminates with matrix material nonlinearity and with initial 
waviness of the fiber was less that the compressive strength of 
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composite laminates with a linear elastic matrix and straight fibers. 
This reference also includes experimental results for compressive 
strength as a function of constituent properties. The authors of 
reference 53 observed changes in the laminate failure mode caused by 
changes in the material properties of the matrix. 

In summary, researchers have studied the failure mechanisms 
for compression-loaded composite laminates for the past twenty-five 
years. " Most of these studies have focused on the short-wavelength 
buckling of unidirectional laminates and on the interlaminar shear 
failures due to lamina imperfections in these laminates. Researchers 
have obtained reasonable correlation between theory and experiment 
using geometrically and/or materially nonlinear analyses. The short- 
wavelength buckling and shear failures of multi-directional laminates 
in compression have received very limited attention. These laminates 
were modeled as the outer-most lamina supported by an elastic 
foundation, and the short-wavelength buckling of the interior laminae 
was neglected. The outer-most laminae always had fibers parallel to 
the loading direction. No general theory exists for analyzing short- 
wavelength buckling and shear failures in compression-loaded symmetric 


multi-directional laminates 
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1.2 Objective and Scope 


The objective of this investigation is to study the short- 
wavelength buckling and the interlaminar shear failure due to short- 
wavelength imperfections of multi-directional laminates loaded in 
uniaxial compression. The inplane shear failure within anisotropic 
laminae is also addressed for these multi-directional laminates. A 
model that focuses on these specific phenomena is presented. The 
model is used to analyze laminates with orthotropic or anisotropic 
laminae. The stacking sequence of each laminate in this study is 
symmetric with respect to the laminate's middle surface. A linear 
analysis is used to determine the laminate stress, strain, and mode 
shape when short-wavelength buckling occurs. A nonlinear analysis 
that assumes initially imperfect (i.e. , initially wavy) laminae is 
used to calculate the laminate stresses and interlaminar strains. The 
initial imperfection for each lamina has the same shape as the 
laminate's short-wavelength buckling mode. A result of this 
investigation is a failure criterion based on short-wavelength 
buckling of the laminae and shear failures within the laminate. 



Chapter 2 


ANALYSIS. 


This chapter presents a model for the short-wavelength 
buckling and shear response of compression-loaded composite laminates. 
The model considers the interlaminar shearing due to short-wavelength 
initial imperfections and the inplane shearing in these laminates. 

The governing equations of the model are derived from first principles 
using an energy formulation. The loading and boundary conditions of 
the model are discussed. The solution of the governing equations is 
described. A linear analysis is performed to obtain the eigenvalues 
and eigenvectors associated with short-wavelength buckling in the 
laminate. A nonlinear analysis is performed to obtain the stress 
distribution in a laminate with initially imperfect laminae. 

2,1 Derivation of the Governing Equations 

2.1.1 Model Description 

The geometry of a typical laminate is shown in figure 2.1. 

The x-y-z coordinate system is referred to as the laminate coordinate 
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system. The laminate is rectangular having length a and width b. 

O 

A +9 fiber orientation in the lamina is indicated in the figure. 

The laminate is symmetric with respect to the middle surface, the z=0 
plane. 

The fundamental element for the model in this study is a 
idealized lamina. The cross-sections of a typical lamina and of an 
idealized lamina are illustrated in figure 2.2. The typical lamina 
has thickness t, and the fibers are modeled as a plate (hereafter 
referred as the "fiber-plate") and the matrix is modeled as an elastic 
foundation (hereafter referred to as the "matrix-foundation"). The 
fiber-plate and matrix-foundation have uniform thicknesses t^. and 

t^, respectively, where 


V 


t 


m 


£(1-V f )t 


( 2.1 ) 


V^. is the fiber volume fraction of the laminate. Combining equations 

( 2 . 1 ), 


t = t + 2t 
f m 


( 2 . 2 ) 
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Idealized laminae are assembled to form the model for a laminate. The 
model for a typical laminate is shown in figure 2.3. The linear 
elastic properties used in this analysis are presented in Table 2.1. 
The fiber-plate properties are typical properties for a graphite-epoxy 
lamina. The matrix-foundation properties are neat resin properties 
from reference 54. The typical lamina thickness is 0.0052 inches. 

Initial imperfections in the fiber-plate are also included in 
the model. Previous authors have referred to these initial 
imperfections as initial waviness of the fibers or as fiber curvature. 
The initial imperfections in this study have the same shape as the 
short-wavelength buckling mode shape for the laminate. 

2.1.2 Fiber-Plate Contributions 


A geometrically nonlinear plate theory is used for the fiber- 

plates. The ith fiber-plate has displacements u^, v and w^ 
in the x, y, and z directions, respectively. The theory for the 
fiber-plates is derived using the Kirchhoff assumptions: 

1 . Normals to the undeformed middle surface of each fiber- 

plate remain straight, normal, and inextensional during 

deformation so that the transverse normal strains e and 

z 
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transverse shearing strains Y and Y in the fiber-plate 

y z xz 

may be neglected; and 

2. Transverse normal stresses a in the fiber-plate are 

z 

small compared to the other normal stress components and may 
be neglected in the stress-strain relations. 

The first assumption leads to the following expressions for the fiber- 
plate displacements, 


u^(x,y,z) = u (i) (x, y) - zw^ (x,y), x 

O 0 

v^(x,y,z) = v^ (x,y) - zw^ (x,y), y 

O 

w (i) (x,y,z) = w (i) (x, y) 


(2.3) 


where u^ , v^ , and w^ are displacements in the x, y, and z 
directions, respectively, for the middle surface of the ith fiber- 
plate, z is the through-the-thickness coordinate of the fiber-plate 
^f — kf 

such that - — S z < 2 * , and subscripts x and y preceeded by 

commas denote partial differentiation with respect to x or y. The 

0 


displacement is the sum of the displacements due to the initial 
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imperfection of the fiber-plate and the displacement due to 

loading ,i.e., 


O 

w (l) (x,y) = (x,y) + w{ x) (x,y) 


( 2 . 4 ) 


The expressions for the normal and shear strains in the ith 
fiber-plate are 


£ 


e 


Y 


(i) 
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0 

X 

(i) 
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♦ lK (1) 
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0 

y 

(i) 
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xy 

xy 

xy 


( 2 . 5 ) 


( ' ) 0 ( . ) 0 ° 

where e 1 , e 1 , and Y ^ are the normal and shearing strains 

* y xy 


for the middle surface of the ith fiber-plate, and k^, and 

x y 


k are the curvatures of the ith fiber-plate. The middle surface 
xy 

strains are defined using the von Karman nonlinear strain-displacement 
relations and include the initial imperfections of the fiber-plate. 
These strain-displacement relations assume that the strains are small 
compared to unity, that rotations relative to the x and y 
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directions are moderately small, and that rotations relative to the z 
direction are negligibly small. The strain-displacement relations for 
the ith fiber-plate are 
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( 2 . 6 ) 


The normal and shearing strains are zero when the fiber-plate is 
unloaded. The curvatures are defined by 


K 


K 


K 


(i) 
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^ XX 
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(i) 

= -w{J } 

y 

i yy 

(i) 

= -2WJ, 1 ) 

1 xy 

xy 


( 2 . 7 ) 


The initial imperfection does not appear in the definition of the 
curvatures since the curvatures are determined by the change in slope 
of the middle surface from the initial middle surface geometry. 
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Each fiber-plate is homogeneous and has linear elastic 
material behavior. A fiber-plate has specially orthotropic material 
symmetry with respect to a principal material coordinate system. This 
coordinate system has axes that are parallel (1-axis) and 
perpendicular (2-axis) to the fiber orientation in the fiber-plate, 

and these axes are illustrated in figure 2.4 for a +9 fiber 
0 0 

orientation (-90 5 9 Si +90 ). The laminate coordinate axes x and y 
are also illustrated in the figure. The constitutive relation for a 
specially orthotropic fiber-plate is 



( 2 . 8 ) 


where {a} and { e } are the inplane stresses and strain, respectively, 
and [Q] are the reduced stiffnesses for the fiber-plate. The reduced 
stiffnesses are calculated using material properties and 
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(2.9) 
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where 


E-j 1 - Young's modulus in the direction parallel to the 
fibers 

E 22 = Young's modulus in the direction perpendicular to 
the fibers 

G 12 = shear modulus in the 1-2 plane 
v^ 2 = major Poisson's ratio 
= minor Poisson's ratio 


The ith fiber-plate has generally orthotropic material 
symmetry with respect to the laminate coordinate system. The 
constitutive relation for the ith fiber-plate is 
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( 2 . 10 ) 


where {a}^^ are the fiber-plate stresses in the laminate coordinate 


system, {e}^ are defined by equations (2.5) and (2.6), and [Q]^ 
are the transformed reduced stiffness of the ith fiber-plate. The 
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transformed reduced stiffnesses are calculated using equations (2.9) 
and 

Q}^ = Q 11 cos J4 9 i + 2(Q 12 +2Qgg )cos 2 e i sin 2 e . + Q^sin^ 

= (Q^ 1 + Q 22 -2 *Qgg)sin 2 0 i cos 2 0 i + Q^Csin^ + cos^) 

= Q.^sin^e + 2 ( Q +2Q, , )sin 2 0 . cos 2 0 . + Q-„cos^9. 

22 ii l 1 2 oo l l 22 l 

®1 6 } = (Q 1 1 ~ Q 12 _2Q 66 )sin 9 i cos3 9i + (Q 12 -Q 22+ 2Q 66 )sin 3 0.cos 0 . 
®26 } = (Q ir Q 1 2' 2Q 66 )3in30 i COS 9 i + (Q i 2 _Q 22 +2Q 66 )sin 
*66 = (Q l1 +Q 22' 2Q 12' Q 66 )sin2e i COs2e i + W 31 ^! + 

( 2.11 ) 

where 0 ^ is the fiber orientation of the ith fiber-plate. For 

O O 

convenience, fiber-plates with e i = 0 or with 0 i = 90 are 

subsequently referred to as orthotropic fiber-plates. Fiber-plates 
with any other 0^^ are subsequently referred to a3 anisotropic fiber- 

plates. 

The potential energy of the ith fiber-plate is the sum of the 
potential energy of the applied loads and the strain energy of the 
fiber-plate. The potential energy of the applied loads is zero for 
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this model as discussed in section 2.1.5 Model Loading and Boundary- 
Conditions. The strain energy of the ith fiber-plate is 
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Using the Kirchhoff assumptions, equation (2.12) is simplified to 
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Using equations (2.5) and (2.10), equation (2.13) becomes 
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Equation (2.14) is integrated with respect to the through-the- 
thickness coordinate z to obtain 
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where A^ and D ^ are the extensional and bending stiffnesses, 
Jk jk 


respectively, of the ith fiber-plate and are defined by 
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The potential energy contribution of each fiber-plate to the total 
potential energy of the model is given by equation (2.15). 

2.1.3 Matrix-Foundation Contributions 


The matrix-foundations elastically support the fiber-plates. 
The ith matrix-foundation is located between the ith and the (i+1)th 

fiber-plates and has displacements u^, and in the x, 

m m m 

y, and z directions, respectively. The displacements of the ith 
matrix-foundation in this model are a linear function of the 
displacements of the adjacent fiber-plates and are defined by 
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where the fiber-plate displacements are defined by equations (2.3), 
and 2 is the through-the-thickness coordinate of the matrix- 
foundation such that -t £ 2 £ t^. Evaluating the fiber-plate 

displacements, equations (2.17) become 
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The initial imperfections of the adjacent fiber-plates do not appear 
in equations (2.18). Deformation in the matrix-foundation is the 
result of displacements due to loading of the adjacent fiber-plates. 

The matrix-foundations in this model have extensional 
stiffness in the z direction and shearing stiffness in the y-z and 
x-z planes. The strain-displacement relations for the ith matrix- 
foundation are 
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Equations (2.18) are substituted into equations (2.19) to obtain 
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Equations (2.20) show that the strains in the matrix-foundation are 
functions of the w displacements of the adjacent fiber-plates. The 
constitutive relations for the ith matrix-foundation are 
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T mxz m mxz 


( 2.21 ) 


where the matrix-foundation is a linear, elastic, homogeneous material 
and 


E - Young's modulus of the matrix 
m 
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G_ = shear modulus of the matrix 
m 

v = Poisson's ratio of the matrix 
m 


The strain energy of the ith matrix-foundation is 



Equation (2.23) is integrated with respect to the through-the- 
thickness coordinate 2 to obtain 
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,( i) 


m 


a b 


n n 2t (1-v ) 
0 0 m m 


ix ' 1 *”**; 1 ’] 2 


at t * ^a )( w ( ;* ,)2 * „<;> 2 + „'f> 2 ) 

4t 2 fm 3 1 y 1'y I'x 1 ’ x 


4t 




(2.24) 


The potential energy contribution of the ith matrix-foundation to the 

total potential energy of the model is given by equation (2.24) where 

i = 1, 2, . .., (N-1 ) where N is the number of laminae in the 

laminate. The two outer-most matrix-foundation layers (figure 2.3) do 

not contribute to the potential energy of the model. The strains in 

the matrix-foundation result from differences in the w displacements 

of the fiber-plates. The outer-most matrix-foundation layers have a 

free surface on one side, and e , Y , and Y are zero for these 

mz rayz mxz 


layers. 


2.1.4 Governing Equations for the Model 


The nonlinear equilibrium equations for the model are derived 
using the principle of stationary potential energy, i.e., the loaded 
model is in equilibrium if its total potential energy n is 
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stationary. This principle is stated mathematically in terms of the 
first variation of the total potential energy as 


6 II = 0 


(2.25) 


The total potential energy of the model is the sum of the 
fiber-plate and the matrix-foundation contributions, or 


N / . \ N-1 , . . 

n = l + l u (l) 

i=i fp i-i m 


(2.26) 


where N is the number of laminae in the laminate. Using equations 
(2.6), (2.7), (2.15), (2.24), and (2.26), equation (2.25) becomes 


[ { l [-(n 4 )+ N (l ! )6u (l) - (N^ )+ N (l ! )6v (l) ] 

I t =1 x x xy’y yy xy x J 

o o 1=1 


+ 2(D ( 1 ) + 2D^ ) )w ( ! ) + 

L 11 1'xxxx 12 66 1’xxyy 22 1 yyyy 

4D ( ! ; ) w ( ! ) + 4D^ ) w{! ) - [N (1 w^)], - 

16 I’xxxy 26 1 xyyy x lx 0 x ’x 

[N (1) (w ( ! } + viV)l, - [N (1) (w ( | ) + w^, } )], - 

y 1y °y y xy lx °x y 


CN (1) (w ( ] ) + w^! } )], + 

xy 1 y °’y x 


2t m (1-v„) 
m m 
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^ [(-§ + 2t t + + w ( ! } ) + 

24 tj,, 2 fm 3 ryy 1 xx 

( "f * 2t f t . * -f Hu 1 ( -yy* »rxx ,] ] 5 ", 0) * 

L [0 ” 1, “'- 1 - 2<D « ♦ 2D «XJU- ^’v ym 

* ^’“rxyyy- ““S'”* "S^x * 

[ f’ ( “ry* ““ y ». y - «<“S”* -‘tyj.,- 

- “x”<“l'” + * 


2t (1-v 2 ) 
m m 


(2w, U) - w' 1H) - „< 1+,) ) 


G t 2 8t 2 

~ C(-| + 2t t + -~)(2w ( J ) + 2w. ( ^ } ) + 
2 f m 3 1 ’yy 1 ’xx 


2 2 

( !i . 2tfV .<£>* »<;;;>* 

«“« »] 5 », U) * W-PL'* MV' 

* »2S W 

[N x" )( “ri* <£>»•, - [»y N><u ry* u '-y )] 'y 


- CN (H) (w'?U 

xy 1 y 0 y ’x 
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xy lx °’xy 




2 v 1 1 

2t (1-v ) ' 1 

m m 


P *.2 o <-2 

u m r , f _ m , , (N) (N) , 

St* C< “ * 2t fV *3* W 1 ’yy Vxx> * 

m “ 


t Hl 

t f . ol . , A m w (N-1 ) (N-1).,i . (N) i , , 

( ~ * 2t fV *3* W 1 ’yy * Vxx Su 1 1 dx * 


N 

l f 
1-1 


N (l) 6u (l) 1 dy + 

N (i ^6u (i) ° 

b 

dx + 

' N (1> <5V H) ° 

x 'o 

xy 

0 

y 


J xy 


> 6 v (i) 


dy - 


M (l) 6w (x) 
x 1 x 


dy - 


M (l) 6w ( , x) 

y i y 


dx 


+ 2M (l) Sw (l) 
xy 1 


] + [M^! J + 2 M (1 ! + N ( 1 ) <wj! } + w<!>) 

J j 1 xx xy y x 1’x 0 x 


G ti 8t 2 

+ N ( W , + w { Q } ) + r-2 [(-§ + 2t t + “T 2 )w 1 ( 1 ) + 
xy 1 y 0 y 4t 2 f m 3 l x 


4* 2 w 


[m ( ! 2M (1 ! + 
y y xy x 


n (,) (. 1 ( : ) * »<:>) * n !,) <w<i>. »!!>> ♦ 

y 1 y 0 y xy l x 0 x 


a » r £ , tt 8t l ( 1 ) 

— [< “ * 2t fV T )w i'y 

m 


4t , * x b N-1 r , . . 

♦ 2tfV >| dx * m j [«';>. 


b 

dx 

0 



a 


2M (l ! + N (l) (w^ )+ w^ } ) +N (i) (w4 )+ w^) + 
xy y x 1 ’x 0 x xy 1 ’y °’y 


2 2 
G t* 8t , . N 

-J2 r (-£ + 2t t + — 2 ) 2 w (l) + 
4t U 2 dZ fm 3 ' 1’x 

m 


4t 


<1 * 2t fV - f K “l < ‘x' ) * "l^; 1 >) ]' SW 1 C1) * * j [”vi 


i(i) 


y y 


+ N (l) (w^ )+ w^) + N (l) (w^ )+ 
xy x y I’y °’y xyl’x °’x 


2 2 
G> t£ 8t ... 

— EJ r (_X + 2t t + — ~)2w 1 + 

4t U 2 f m 3 1’y 

m 


4t 




[(,«>* 2M (N > * » (N) (W«>* v'?>) * 
J 1 xx xy’y x lx 0 ’ ~ 


G t 2 

m r / r 


8t c 


N (N) (w (N) + w (N) ) + —2 [(-£ + 2t t + -^)w (N) + 
xy 1 ry °’y ; 4t m U 2 m 3 ; I’x 


2 2 

(!l , 2, t * ^J5)„ <N -' h s „«> \y 

1 2 f m 3 1 ’X J6W 1 _ y 


[M ( ^+ 2M (N ) 

L y y xy x 


♦ * W dV+ w?) + «iV) + 

y 1 y ° y xy l x 0 x 


r 4-2 

_jn r r— + 2t t + ~2)w (N) + 
4T C( “ + 2t f t m ~ )w l’y 
m 



fir ♦ 2t t * J)„< N - , h s „ <,<> b d X 
1 2 m 3 ; Vy J 1 0 


(2.27) 


where the stress resultants for the ith fiber-plate N^, and 

x y 

are defined by 
xy 


N (i) = ASj ) e (i) ° + A$J ) e (1)# + A 1 ( j ) Y (i) ° 
x 1 1 x 1 2 y 1 6 xy 


N (l) = A ( J } e (l) + A^ } e (l) + A^\ (l) 
y 12 x 22 y 26 xy 


( 2 . 28 ) 


N (l) = A ( ^ } e (l) + A^ } e (l) + A^\ (l) 
xy 16 26 y 66 xy 


and the moment resultants for the ith fiber-plate M ^ , and 

r x y 


are defined by 
xy 


-(D (l) w ( ^ + + 2D ( ^ } w ( ^ ) 

11 1 ’ xx 12 1 ’ yy 16 1 xy 


12 1’xx 22 1’yy 26 1’xy 


(2.29) 


16 1 xx 26 1 yy 66 1 ’xy 


The summations and integrations in equations (2.27) can be 
interchanged since the series of potential energy contributions 
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consists of continuous functions and is uniformly convergent over the 
region OSxSa, 0 i y To (see reference 55, p.589). 

Equation (2.27) is satisfied when the integrand of the area 
integral and when the integrands of the boundary condition integrals 
equal zero independently. The governing nonlinear partial 
differential equations are obtained from the area integral and are 
stated below: 


for i = 1 , 2, . . . , N 


N ( ^ + N (l J 

x x xy y 


= 0 


+ N (l) , 
y y xy x 


( 2 . 30 ) 


for i = 1 


d (1) k (1) ] - n ( 1 W!> ♦ w<!> ) - 

L 1 J X 1 XX 0 XX 


N^^w^P + ■,<]> ) - 2N (1) (w ( ] ) + w^ 1 ,' 1 ) 
y T yy 0 yy xy 1 xy 0 xy 

+ £*(w (1) - W (2) ) - G*(w (1) + w (1) ) - 
a ^ W 1 W 1 } W 1 yy Vxx ; 


G*(w ( f 5 + w (2) ) = 0 
2 r yy 1 xx 


(2.31) 



where 
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for i = 2, 3 N-1 

D 


(1) [w, (1) ] - n (1) (w<; ) ♦ > - 


1 ' XX 


, (1) 

O’XX' 


N (l) (w^ + ) - 2N (l) (wJ^ + ) + 

y T yy 0 yy xy r xy 0 xy 

E*(2w (i) - w (1 ' 1) - w (i+1) ) - G*(2w (i) + 2w (1) ) 

t w 1 w 1 ; V Vyy Vxx ; 


G*(w4~ 1)+ w4~ 1)+ w4 + 1)+ w4 + 1 - ) ) = 0 

2 1 ’ yy 1 xx 1 yy 1 ’ xx 


(2.32) 


for i = N 


^N)[ <M)] . N <») ( (N) , „(») , . 

L 1 J X 1 XX 0 ’ XX 


» (N) („«) * „«> ) * „«> ) . 

y 1 yy 0 yy xy Txy °’xy 


* (N) (N-1). *, (N) ^ (N) . 

E (w, - w. ) - G„ (w, , + w. , ) 

i i i i yy i xx 


G*(w<“> . «<?> ) - 0 
2 1 ’ yy 1 ’ xx 


(2.33) 


D^[ • ] is an operator defined by 


D (1) [w] = D.\ i; w, + 2(D 1 V ^+ 2D^t J )w, h 

L J 11 ’ xxxx 12 66 xxyy 


.(i) 


(i). on(i) 


D^w, + 4D^ P w, + 4D^w, (2.34) 

22 yyyy 16 xxxy 26 xyyy 
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and the constants are defined by 


# 

E 


E 

m 


2t (1-v„) 
m m 


ir<r * 2t rV 

m 


8t 


2 2 

G m *7 4tf 

_J2(— £ + 2t t + — - — ) 
4t ^2 m 3 

m 


(2.35) 


The equations for the boundary condition are obtained from the 
boundary condition integrals and are stated below: 

for i = 1 , 2, . . . , N 


N (1) 5 U (1) 

X 


N (1, SV (I) 

y 


= o 


o 

0 ■ b 


= 0 


m<V‘’ 

X l x 


N (l) 6u (l) 

xy 


N (l) 5v (l) 

xy 


= o 


M (i) 5w 1 ( f ) 
y i y 


o 

0 ■ a 
0 
b 


= 0 


= 0 


2M (l) 5w^ l) 

xy 1 


(2.36) 
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for i = 1 


[M (1 , *+ 2M (1 ! + w<!>) + 

1 x’x xy y x 1'x °x 


N 


( 1 ),..( 1 )| 

xy 


(2)i,„(1) 


(w u; + w u; ) + G w u; + g w l liw 
V Vy °’y J °1 w 1’x Wx Jd 1 


= 0 


[M ( l } + 2M (1 ! + N (1) (w ( ] } + w[! 5 ).+ 
1 y y xy x y I’y °’y 


N (1) (w (1 /+ w<!>) + G*wj] J + G*w ( ^]6w (1) 
xy 1 x 0 x 1 1 y 2 1 y J 1 


for i = 2, 3 N-1 


[m ( ^ + 2M (i ] + N (l) (w ( * } + w^ } ) +N (l) (wj/ ) + 

L xx xy’y x I’x °’x xy 1 y °y 




= 0 


0 


[« ( ;4 2h ( 1 ! * w<») ♦ n (1, (w!;4 »!;>> 

L yy xy’x y I’y °’y xy 1 x °x 


* 2C *“ry * “ry ) )I<"i (l> 


= 0 


for i = N 


[«<»>.. 2M ( "! . N (N) (w'“4 „<?>) + „«>) 

L x ’ x xy’y x I’x 0 ’ x xyl’y °’y 


+ gV N) + G*w (N ’ 1) l6w (N) 

hVx U 2 W 1’x J6W 1 


= 0 


(2.37) 


( 2 . 38 ) 
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[M ( ? } + 2M (N ! + w<?>) + w<?> 

1 y y xy'x y 1 * y 0 * y xy l’x °’x 

* (N) * (N-1)i.C N )| b . 

G 1 W 1 ’ y + G 2 W 1'y ^ 6W 1 L " ° 


(2.39) 


2.1.5 Model Loading, Boundary Conditions, and Solution 
Procedure 

The model is loaded in compression by uniform end shortening. 
No external loads are applied to the model (hence, the potential 
energy of the applied loads is zero). This loading enters the 
analysis as a u displacement boundary condition for each fiber-plate 
and is expressed as 


u (l) (0,y) = 0 

0 

u (l) (a,y) - -u 0 ,1-1, 2, .... N (2.40) 

The u and v displacement boundary conditions are chosen to 
assure that the model remains rectangular during loading. A 
rectangular shape is typical for composite laminates in structural 
components. These boundary conditions for each fiber-plate are 


expressed as 
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The v displacement boundary conditions are the result of Poisson 
expansion during loading. These boundary conditions prevent an 
anisotropic lamina from deforming into a parallelogram-type shape. 
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The model is simply supported on the unloaded edges. These 
simple-support boundary conditions are expressed as 

at y = 0 and y = b 

Wl (l) = 0 1 } = 0 ,1-1,2, ...? N (2.43) 

The model with initial imperfections is loaded in the x 
direction and deforms into a short-wavelength buckling mode. The 
half-wavelength of the model's mode shape X is a laminate property, 
and the calculation of X is discussed in section 2.2 Linear 
Analysis. This half-wavelength is assumed to be much smaller than the 
laminate's length or the laminate's width. The fiber-plates in the 
model are assumed to behave as elastically supported semi-infinite 
strips causing the w displacement and moment boundary conditions at 
x=0 and at x=a to have a negligible effect on the short-wavelength 
buckling behavior of the model. The w displacement and moment 
boundary conditions for this model are expressed along node lines of 
the mode shape and are 
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at x » 0 and at x = X 

w{ l} = 0 M^ l} = 0 , i = 1, 2, .... N (2.44) 

where the x-axis is parallel to the x-axis and is illustrated in 
figure 2.5. 

The Kantorovich method [56] is used throughout this study to 
obtain solutions to the governing equations that satisfy the boundary 
conditions. Each unknown is expressed as a kinematically admissible 
series of the x coordinate. An unknown function in this study has 
the general form 

f(x,y) = f 0 (x) + f 1 (y) + f g (y) sin ~ + 

f c (y) cos 22 (2.45) 

where f 0 and f^ are known functions determined from the linear 

prebuckling analysis. Equation (2.45) is substituted into the 
governing equations, (2.30) to (2.33), and boundary conditions, 
(2.41), (2.43), and (2.44), to obtain ordinary differential equations 
for the unknowns f and f . 

3 C 
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2.2 Linear Analysis 


The nonlinear governing equations are linearized to determine 
the laminate end shortening when short-wavelength buckling occurs. 

This end shortening is also referred to as the critical end shortening 
( u °)cr* The linear * zecl equations are also used to determine the 

short-wavelength buckling mode shape. The linearized equations are 
derived using the adjacent-equilibrium criterion (see reference 57, 
p.27). The fiber-plates have no imperfections in the linear analysis 

(i.e. , w[^= 0). 


Adjacent equilibrium configurations for the fiber-plates are 
investigated by adding small increments to the fiber-plate 
displacements. The relationships between the displacements prior to 


k (i) (i) . (i) 

buckling, u , v , and w , 

cl a 3. 

equilibrium configuration, u^\ 


and the displacements of an adjacent 
and w (i) are expressed by 


u (i) =u (i) + eu 
a 


V (1 >. v li) t E 7 
a 


u (i). w (‘). c ; 

1 a 1 


(2.46) 
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where the quantity e is infinitesimally small. The prebuckling w 
displacement and its derivatives are zero for the flat ith 

3 . 

fiber-plate loaded by uniform end shortening. 


Equations (2.46) are substituted into equations (2.30) to 
(2.33), and the terms of like powers of e sum to zero for arbitrary 
e. The zeroth-order terms in e combine to give the linear 
prebuckling equilibrium equations for the ith fiber-plate 


N 

N 


(i) 

+ N 1 ' 1 ] 

X 

o 

X 

xy° 

(i) 

y°’y 

+ N (l ^ 
xy° 


y 


X 


= 0 
= 0 


(2.47) 


where 


N 

N 

N 


(i) 

x° 

(i) 

y° 

(i) 

xy° 


A (l) u4 ) + A ( J ) v4 ) + A 1 4 ) (u4 )+ v ( ^) 

11 ax 12 ay 16 ay ax 

12 a’x 22 a’y 26 ay a’x 


♦ A< 1) v<1 ) * a»)(u ( !>» v<;)) 
16 ax 26 a’y 66 ay a’x 


(2.48) 


The prebuckling displacements that satisfy the governing equations 
(2.47) and the boundary condition equations (2.40) and (2.41) are 
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u ^ i) (x,y) 


- /Xs 


V l l)<X>y) ' T 
A 22 


a t 

(~)(y 


, u 

b) _I 

2 a 


(2.-49) 


The first-order terms in e combine to give the linear postbuckling 
equilibrium equations for the inplane loads in the ith fiber-plate 


N (l ! + N (l) , 
xl x xyl y 


N + N (l) , 
yl y xyl x 


= 0 


(2.50) 


where 


»<?)- A;, 1 );! 1 * . a<‘>7<‘> . 7, (1 >) 

xl 11 x 12 y 16 y x 

(i) (i)-(i) .(i)-(i) . ( i) /~( i), — ( i ) n 

"yl ' *12 u ’x ft 22 v ’y * A 26 <u ’y * v >x > 

(i) (i)-(i) (i)-(i) .(i),-(i). — ( i ) n , p r.v 

N = A,, u, + A„ r v, + A rr (u, + v, ) 12.51.) 

xyl 16 x 26 y 66 ’y x 


and the linear stability equations for the fiber-plates are stated 


below: 
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where 


for i = 1 


q(i ) r— ( 1 ) 


(1) w (1) - N (1) w (1) 
x° 1 ’xx y° 1 ’y 

*-(1) -(2), *,-(1)' -0) 


[w; U ] - - N - 2N (1 ]w 1 ( ! ) 

1 J x° rxx y° 1 yy xy° 1 


y u i 'yy xy° 1’xy 


+ e (*;■'- w ;^ y ) - g.( W ; ' + w; ' ) 
1 i 1 1 yy 1 xx 


G*(w (2) + w (2) ) = 0 
2 1 yy 1 xx 


(2.52) 


for i = 2, 3. ...» N-1 

rr(i)r-(i)i _ „(i)-(i) _ „(i)-(i) 


t < i )— ( i) 


D vi; [w l '] - N o wj , ; - N [ yw y - 2N — w- 
L 1 J x° 1 xx y° 1 yy xy° 1 xy 


+ * (2w l W 1 W 1 } 2G l ( Vyy + Vxx } 


1 


'1 


♦ g!^!" 0 * w;r 1} + w;i +i) + w;i +i) ) = 0 

2 I’yy 1 xx 1 ’ yy 1 xx 


(2.53) 


for i = N 


D (N >R«>) - N ( " ) w« ) - N C :>5<?> - 2N<«W<? ) 

L 1 J x° 1 xx y° 1 yy xy° 1 xy 


* (N) -(N-1), *-(N) , — (N) . 

+ E (w, - w. ) - G, (w. , + w. , ) 

1 1 1 1 yy 1 ’ xx 


21’ yy 1 ’ xx 


(2.54) 


and are determined from equations (2.48) 

x° ’ y° xy° H 


The higher-order terms in e are neglected since e is small 
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2.2.1 Orthotropic Fiber-Plates 

The short-wavelength mode shape for an orthotropic fiber-plate 
is assumed to have half-waves that are normal to the direction of 
loading. These half-waves are referred to as "normal waves" and are 
illustrated in figure 2.5a. The linear ' stability equations and 
boundary conditions for a model composed of orthotropic fiber-plates 
are satisfied by 


:(i) 


— (i) . 

= w sin 
s 


irx 

X 


sin 


iry 

b 


(2.55) 


where w^^ is the amplitude of the normal wave for the ith fiber- 

TTX 

plate. The coefficient of the cos — - term is zero (see equation 

A 

(2.45)), and the extensional stiffnesses, A^ ^ and A^g^, and the 

bending stiffnesses, and are zero f ’ or orthotropic fiber- 

plates. Equation (2.55) is substituted into equations (2.51) to 
(2.53) to form a system of equations for calculating the model’s 
critical end shortening and short-wavelength buckling mode shape. The 
system of equations is of the form 
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['d (N -’ >. 2M 


k 


N-1 



-(N-U 


Nw 


(N-2) 



0 




»s (N - 1) 


- 0 ( 2 . 56 ) 


for arbitrary x and y where the constants are defined by 


n(i) 


>( i ) ^ ™(i) 


" D H T + 2(D 12 + 2D 66 )( A } ( b> + °22 


M = E + C,[(¥ + (?) d ] 
1 b X 


* * 17 ? TT ? 

e - G_[(£r + (f) ] 

2 b X 


k . r A (i) - OS A (i) ](5) 2 
i 1 11 T 12 J V 

A 22 


r A (i) - -IS a < 1 >1(S ) 2 

1 12 „T 22 M b ; 
A 22 


(2.57) 


Non-^trivial solutions to equations (2.56) exist if the determinant of 

the matrix formed by the coefficients of is zero. The smallest 

s 

value of the normalized end shortening u 0 /a for which this 
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determinate is zero is the critical normalized end shortening 
(u 0 /a) cr . The end shortening is normalized by the laminate length, 

a. The value (u 0 /a) cr is an eigenvalue of this system of equations, 

and the corresponding eigenvector is the short-wavelength buckling 
mode shape of the model. The wavelength of the mode shape is 
determined by minimizing (u 0 /a) cr with respect to the wavelength. 

2.2.2 Anisotropic Fiber-Plates 


The short-wavelength mode shape for an anisotropic fiber-plate 
can have half-waves that are oriented at an angle <j> to the 
orientation of a normal wave. These half-waves are referred to as 
"skewed waves" and are illustrated in figure 2.5b. The form of the w 
displacement that satisfies the linear stability equations and 
boundary conditions for a model composed of anisotropic fiber-plates 
is 

— (i) (i), A . ttx t (i), . irx ,- Q . 

w 1 = w s (y) sin ~ + w Q (y) cos — (2.58) 

The eigenvalues and eigenvectors for these fiber-plates are obtained 
using the Galerkin method [57]. The functions of y for the ith 

fiber-plate, w^ and w^ , are expanded in a Fourier sine series, 
s c 
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These series satisfy the boundary conditions for w j ^ at y=0 and at 
y=b and are expressed as 


(i) r — (i) . 
w = ) w sin 
s L , sra 
m=1 


m-rry 

b 


(i) r -U) 
w = ) w cos 
c L , on 
n=1 


niry 

b 


(2.59) 


where M is the number of terms in the series and is an even number. 
Equations (2.58) and (2.59) are substituted into equations (2.51) to 

(2.53), and the series expansions for w^ and are solutions 

S G 

miry 

to these equations if the equations are orthogonal to sin -jp and to 
niT v 

sin * for each m and n, respectively. The orthogonality 
conditions for the ith fiber-plate are stated below: 


for i = 1 


(d 0> [w 5 ,j ] - > - 

J L 1 J x° 1 xx 


dfcO) _ .(IHI) . 2B o>u<p 

y o 1 9 yy xy° 1 xy 


+ E *(w (1) - w (2) ) - G*(w (1) + w (1) ) - 
h ^ W 1 W 1 ; Wyy Vxx' 

,*,—(2) t -(2) 


G 2 (w 1 ’ yy + W 1’xx^ sin ? dy = 0 m = 1, 2 M 
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for i =2, 3> . . . , N-1 


[D j - N x0 w 1>xx - N y0 w Vyy - 2N xy0 w 1>xy 


+ * r? -(D — ( i“1 ) -( i + 1 ) -> _ + -U) v 

+ E (2w i ~ w t w i ) 2G 1 W 1 ’yy W 1 ’ xx 5 

+ g*(w^'!. )+ w^~1 )+ w^lV + sin ^ dy = ° 


2 1 yy 1’xx 1’yy 


m = 1 , 2, . . . , M 


for i = N 


f-Wr-WT n (N)-(N) _ (N)-(N) _ ? „(N)-(N) 

(D [ Wl J " N x0 w 1>xx - N y0 w 1>yy 2N xy0 w 1>xy 


m = 1 , 2, . . . , M 


for i = 1 


x° W 1’xx y° W 1’yy xy° 1 xy 
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“! ? xx» •*» ¥ ^ 


1' Tyy 1'xx‘ 

0 n = 1,2 M 


for i = 2, 3, .... N-1 
b 


- 2N^v/^ 

1 x° 1 ’ xx y° W 1’yy xy° w 1’xy 


+ E*(2w 1 (i) - wj 1 " 1 ’- w (i + 1) 


;(i) ^ — ( i) 


1 


) - 2G.(w + w r ) 

1 1 yy I’xx 

r */~ (1-1). — ( i— 1 ) -(i+1) — (I+1)>'i . mry , 
G 2 W 1 ’yy * Vxx * W 1 *yy * Vxx >) 310 b ay 

n = 1 , 2, . . . , M 


for i = N 


(D«> - N ( »W> - 2N (N ^«) 

1 J x° 1 xx y° 1 yy xy° 1 xy 


1 


* £*(='»)- If-')) -g*(5<?> * ;<« ) - 

1 1 1 1 yy 1 xx 


*-(N-1) . -(N-1 ) x i . niry . 
G 2 ( Vyy + W 1'xx )] Sin b dy = ° 


n = 1 , 2 , ...» M 


( 2 . 61 ) 


Equations (2.60) and (2.61) are integrated to obtain 2*M 
simultaneous homogeneous equations for each fiber-plate in terms of 
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the 2*M undetermined coefficients and w^. The 

sm on 

simultaneous equations for the ith fiber-plate can be expressed as 


L odd 

H 

0 

0 

(i) 

r > 

^ w s^odd 

(i) 

H 

L 

even 

0 

0 

< 

■ ^ W c^even 

= {0} 

> 

0 

0 

L 

even 

-H 


^ w s^ even 


0 

0 

-H 

L odd 


^ w c^odd 

V. J 

(2.62) 


where the submatrices for the ith fiber-plate are defined by 


* U o 


L r k i ( r } 


(i) 


* U o 


0 L -k ( — ) 0 

3 3 a 


L (i) 

L odd 


* U o 


0 L 5 ' k 5 ( r ) 


(2.63) 


* U o 

L (M-1 ) -k (M-1 ) ( a" } 



* U o 


L (l) 

even 


L 2 _k 2 ( r ) 0 


* U o 


0 L 4 -k 4 ( r ) 0 


* U o 


0 L,-k,( — ) 

o 6 a 


* U ° 
L (M)' k (M) ( a } 


(i) 


H 12 h 12 ( a ) H l4 -h l4 ( a ) 


1 M 1M a 


H (i) = 


H 32~ h 32 ( a } H 34 h 34 ( a ) 


3M 3M a 


H (M-1 )2 _h (M-1 )2 ( a } 


H (M-1 )M h (M-1 )M ( a 5 


(2.64) 


(i) 


(2.65) 


with 
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l(i) = D (l) (J) 4 + 2(D ( J } + 2D^ ) )(2) 2 (2Z) 2 + d^ } ( 2I) 4 
m ii A 1 2 Ob A D 22 b 


+ c t (i) [E* + G*[ (2) 2 + (SI) 2 ]] 


.( i)_ 
m 

c <*“’ 

A 12 

a t 



A 22 

.( i)_ 
inn 

:« d »> 

(I) 3 + 

V 

U) = 

<- 

a t 

12 

mn 

A T ' 



A 22 

,( i)_ 

1 , i = 

1 or 


^2 H P + (A 12 t )( b >' ] 

A 22 


4mn 


26 V\b 


, / 2 2 v 

b(m - n ) 


KZ 'Hr) 


8mn 


. / 2 2 . 
b(m - n ) 


= 2, otherwise 


and the sub-vectors for the ith fiber-plate are defined by 


s3 

(w J w p 

s odd ^ s5 


w 


si 


w 


w 


(i) 


s(M-1 ) 


<*> 

W s2 


s4 


(w ) ^ = J w , y 
s even > s6 r 


w « 
sM 

V. J 


2 . 66 ) 


(2.67) 
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The orthogonality of the trigonometric terms causes coupling between 

the (w ) ^ and the (w coefficients and between the 

s odd c even 

(w and the (w ) ^ coefficients. Equations (2.62) lead to an 

s even c odd 

eigenvalue problem of order 2#M*N for an N-lamina model. This 
2*M*N system has only M*N unique eigenvalues. The unique 
eigenvalues are obtained form equations for the fiber-plates of the 
form 




= { 0 } ( 2 . 68 ) 


The system of equations for the N-lamina model are shown in figure 2.6 
and are labeled equations (2.69). The undefined sub-matrices in 
equations (2.69) are 
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N 1 0 0 

o n 3 0 

0 0 N c 

0 

0 0 0 

N 2 0 0 

0 0 

0 0 N,. 

o 

0 0 0 

with 

N = E* + G*[(«~) 2 + (?) 2 ] (2.72) 

m 2 d A 

The N 3ub'matrix is the result of coupling between adjacent fibers 
plates through the matrix-foundation. Non-trivial solutions to 
equations (2.69) exist if the determinant of the matrix of material 
properties is zero. The smallest value of u 0 /a for which this 

determinate is zero is (u 0 /a) c ^, and this value is an eigenvalue of 
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equations (2.69). The eigenvector corresponding to (u 0 /a) cri is used 

to determine the short-wavelength buckling mode shape of the model. 
Similar to the model with orthotropic fiber-plates, the wavelength of 
the mode shape for a model with anisotropic fiber-plates is determined 
by minimizing (u 0 /a) cr with respect to the wavelength. 

2.3 Nonlinear Analysis 


The governing nonlinear differential equations are used to 
obtain displacements, strains, loads, and stresses for a model with 
initially imperfect fiber-plates. The initial imperfections for all 
fiber-plates are the same and have the shape of the short-wavelength 
buckling mode for the model. The Kantorovich method used in the 
linear analysis also is used to obtain solutions to the nonlinear 
equations. The unknown functions are expanded in kinematically 
admissible trigonometric series in the x-coordinate, and the governing 
equations become nonlinear coupled ordinary differential equations. 

An approximate solution to the nonlinear coupled ordinary 
differential equations is obtained using an algorithm developed by 
Lentini and Pereyra [58]. The algorithm is based on Newton's method 
and can be used to solve a system of simultaneous first-order 
nonlinear ordinary differential equations subject to two-point 
boundary conditions. The algorithm uses finite differences to solve 
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the differential equations, and adaptive mesh spacings are 
automatically produced so that mild boundary layers are detected and 
resolved. The difference between successive approximations for each 
unknown is calculated. The algorithm iterates until this difference 
is less than a specified tolerance. 

2 . 3 » 1 Orthotropic Fiber-Plates 

The nonlinear analysis for models with orthotropic fiber- 
plates is similar to the analysis used by Stein [59] for the 
postbuckling behavior of orthotropic plates. The displacements in the 
present analysis for the ith fiber-plate are expressed as 


(i) 


-u 0 x 


( l) , , . 2irx 

— + u (y) sin -7— 

a s J \ 


(i) (i) , . (i)/\ 2 ttx 

r = v 0 (y) + v (y) cos -r— 

G A 


(2.73) 


(i) /X . TTX ( i ) . v . TTX 

* = W 0 (y) sin — + w s (y) sin j" 


Equations (2.73) are substituted into equations (2.6) and (2.7) to 
obtain these expressions for the strains and curvatures: 
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(i) 


• \ 2 


0 2ir (i) 2 ttx 1,ir.2 (i)*", 2irx s 

- + — u„ cos — r— + t(t) w (1 + cos - 7 -) 


4 V X' 


1 ,ti\2 ( i) . 2 ttx , 

♦ 2 ( »> "."s (1 * 003 T 1 


(i) (i)' ( i) • 2 ttx x 1, ( i) 1 >2, 2itXx 

e y = v 0 + v c cos — - 4 ( w s ) H - cos — ) 


1 ' ( i ) ' , 2 ttx n 

+ - w 0 w g (1 - cos — ) 


v (i) (l)' . 2 ttx 2 it (i) . 2 ttx 

Y = u sin — — - — v sin — — + 
xy s X A c X 


bH l) ^ y 


sin 


2-rrx , 1 , ir 


+ ;(f)w D w 


( i) ’ 


sin 


2ttx 


1 ytr s ' ( i ) . 2 ttx 
* 2 X W ° W S S1 " T 


( i) ,ik 2 ( i) . ttx 

K x ■ h> " s S1 " r 


to 

'7 

(i) 


- w 


(i) 


sin 


TTX 


0 ,ir, ( l) ' TTX 

k ■ ' = -2(-)w cos — 
xy X s X 


(2.74) 


where ( ) denotes differentiation with respect to the y-coordinate. 
The stress and moment resultants are expressed as 
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where 


M 


M 


[ (i) = 

X 

N^ } (y) 

+ »re <*> 

2irx 

COS — “ 
A 

(i)_ 

y 

Ny 0 i} (y) 

+ <o )( ^ 

2 ttx 

cos 

( (1) = 
xy 

N (l) (y) 
xys J 

2irx 
sin — 

A 


i (i) = 
X 

M (l) (y) 
xs J 

irx 

sin — 

A 

- 

I (i) = 

y 

M (l) (y) 
ys J 

TTX 

sin r— 
A 


, (i) = 
xy 

M (l) (y) 

xyc 

TTX 

COS r— 
A 



,U) .(i), 


u. 


,x 2 


n ;:'- r f i'i'M 11 * * 


V U ') 2 * 1 w'» (1), l 

A 1 2 L 0 it 1 s 1 2 ° 3 J 


4 U * ;<i>M 1>2 * l<?> 2 “.““^ * 


12 L c ir s ' 2 ° s J 


(i) = A (i) Cv (i)' + l (w (i)’)2 + i w ' w (i) ? 3 + 

V A 22 LV ° IPs ; 2 W ° s J 




N (1) . A (1) rv (i) ’- i(w (1) ') 2 - 1 wV 1 *'] . 
ye 22 L o 4 l a ' 2 » s J 


A i ( 2 >[ r “a 11 * 5 ( j ) 2 “s 1>2 * l ( ! )2 "."“ !] 


(2.75) 
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» (i) - fl v (1) . 1(1) vW 1 ’’. 

xys 66 5 X c 2Xss 


1 /if 


( 7 ) w n w 


(i) ' 1/1 


(i) 




2 X ' 0 s 

M (i) = D (i) ( I)2 w (i)_ D (i) w (i)" 
xs 11 X s 12 s 


M y“" - D 22 ) “i 1) "" “I"'!''""’ 


M (l) = -2D^(;f) w (l) ' 
xyc 06 X s 


(2.76) 


Equations (2.74) to (2.76) are used in equation (2.25) (principal of 
stationary potential energy), and the resulting equation is integrated 

over 0 < x < X. For arbitrary 6 u^\ 5v^\ 6 v^\ and 6 w^\ the 

cL d- S 

principle of stationary potential energy requires that the following 
differential equations be satisfied: 


for i = 1 , 2, . . . , N 


fi *u>- 0 

X xc xys 

N ( o ),= 0 
y° 


flN (i) + N (i) ’=0 
X xys yc 


(2.77) 
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for i = 1 


5 ( ! )N xyi <sl1) * u » ) ‘ G *C-8 <n '* <j) 2 "s ,>] ' 


AS S S 


( 2 . 78 ) 


for i = 2, 3, ♦ ♦ . ,N 

H>0 6<1) * ■ 2 <[-8 (1> '* ( !> 2 "s 1>] - 

2 AS AS 


- E*(2w (l) - w (l 1 w (l + 1 ) ) 
s s s 


( 2 . 79 ) 


for i = N 


*r- ».>- 

2 ( X )M xy3 <6<N) * “(> ' 


G*:-e (N - , )% - eV" 1 - (2.80) 

2 A S S S 


where the previously undefined functions for the ith fiber-plate are 
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V (i) = -M (l) + ~ M (i) - (N^- \ N (l) )(g (l) + w ;> - 
y ys X xyc y° 2 yc 0 

l(7)N (l) (w (i) + w ) 

2 X xys s 0 


( 2 . 81 ) 


g (l) = w (i ^' 
s 


The boundary conditions for this energy formulation are 


N (i) 6u (i) 
xys s 


N (l) Sv (i) 
yc o 


= 0 


= 0 


yO o 


M ( x) 6w ( ' 
ys s 


= 0 


= 0 


( 2 . 82 ) 


V (l) 5w (l) 
y s 


= 0 


i — 1 > 2 , ...» N 


The model remains rectangular during loading and is simply-supported 
on the unloaded edges. These boundary conditions are expressed as 


at y = 0 and at y = b 

.U) „U) „U>_ „(D_ M (i) 


u 


v = N „ = w v = M' = 0 
o y° s ys 


(2.83) 


The equation (2.77) and boundary condition equation (2.83) for N 


(i) 


are satisfied by N^ = 0 for 1 i i S N throughout the region. 
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The first-order governing equations for a model with 
orthotropic fiber-plates can be expressed as a function of eight 

unknowns for each fiber-plate, and these unknowns are u^*\ 

s c 

w (i) N^, and The other unknowns are 

s yc xys ys y 

expressed as functions of these eight. The eight differential 

equations for each fiber-plate are obtained by re-organizing equations 

(2.74) and (2.77) to (2.81) and are stated below: 


for i = 1 , 2, . . . , n 
A 66 


,(i) ' 


N 


(i) 


(i) 


c_ 1 fl (i) 2 + 1 » (i)_ n (i) + 

try 4 8 + 2 w ° 3 t n [ r u s + 

22 22 


1 . -JT. 2 (i) 2 Un.2 (i), 
4 X W s + 2 X W = W s ] 


M (i) n (i) 

' " D U) D U>V W S 

22 22 


e (i) 

s 


(2.84) 
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.Cl) 




N (i) '= (i) + ( A (i).il2 )[ | I (i) + 


xys 'X '\(i) yc 
k 22 


‘11 (i) ' U X s 

22 




N (i) ’= -( — ) N (l) 
yc X xys 


« (i) ‘ 

ys 


-v (l) - 4D^ } (r) 2 s (l) + i w ’ a ) ~ 


'66 K X 


2 yc 


5 N (l) (7)(w (l) + w ) 
2 xys X s 0 


for i = V 


.(D 


.0 )‘ 


v y 1), ‘ ' ( j >2 i“jfr H ys > * (D , ( ; ) - 7 tt )( x )2 "1 1>1 - 
U 22 22 


A 


(1) 


22 


, (A n>_ iilL):- 2 2l * Si u <d* 

U 11 .(1) ;L a X s 

A 22 




M U) D (1) „<2) 

ys ( 12 ,,TT,2 (1)-. r *r ys. 

G i c rtTT (1 * jiT u x ) “s 3 G 2 [ n ( 2 ) 

u 22 22 22 
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D 


(1 - 


(2) 

12 wir.2 (2)-, 

rci)^ w s j 
22 


(2.85) 


for i =2, 3, ...» N-1 

v (i) '= -(t) 2 [-4tt M ( i )+ (dJi 0 - ^-)(7) 2w 1 U ] - 


‘A 7 L n (i) ' ys 
J 22 


'11 D (i) 
22 


A 


(i) 


A 22 


A (i) 2U 

* < u - 4fr>t- -r * f 1 

A 22 


wf>], -eW‘>- .“-’l- 


s 
(i) 


„ M 
*- ys 


1 D (i) 
22 


(i) 

(i) ; V *s 

22 


M 


s 

(i-1) 


- 2G. [-4rT + (i - -nr)(r) t w' i/ ] ~ g*[-^ 


2 D (i_1) 
22 


D ( i-1) M (i + D 

(1 - Ji )( I } 2 (i-i ) + .ys..,. + 

U (i-l)'V W s (i+1) 

22 22 


D 


(1 


(i+1 ) 

12 XriV 1 * 0 : 


D 


(i + 1)'V s 
22 


( 2 . 86 ) 
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for i = N 


>(N) 


,(N)‘ 


(D (N) - -!i-)(I) 2 v / N >l - 
y V l n (N) ys ^ U 11 (N) ; V s J 

22 22 


IC^ 00 * 

A 22 


A (N) 2U 

+ (a (N). Ji_ )C _ _I + 2w (N) 
^11 . (N) >l a X U s 

22 


3,^2, (N) 2,-,, * (N) (N-1), 

(w g + w 0 )]} - E (w g - w g ) - 


M 


(N) 


,(N) 


(N-1 ) 


*r .JfsL + (i - . 1 2 „.) (H) 2 y/ N ) ] - g*[-XS + 


1 V N) 

22 


( N ) X ' "s 
22 


2 L (N-1 ) 
22 


D 


(N-1 ) 


<’ -*T7>W°> 

°22 


( 2 . 87 ) 


The function of y for the initial imperfection w is expressed as 


w 0 (y) = w sin 


n 


( 2 . 88 ) 


where w is a constant that is input to the analysis. 
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2.3.2 Anisotropic Fiber-Plates 


The nonlinear analysis for compression-loaded models with 
anisotropic fiber-plates is similar to the analysis used by Stein [60] 
for the postbuckling behavior of orthotropic plates loaded in combined 
compression and shear. The equations in the present analysis for the 
ith anisotropic fiber-plate are general forms of the equations for the 
ith orthotropic fiber-plate. The displacements for the ith 
anisotropic fiber-plate are expressed as 


l(l) = *” + u o (i) (y) + u^ i} (y) sin + 


(i), 2-rrx 

u q (y) cos 




2tt x 


(y) + ; (y) sin *2. 

S A 


+ ^ l} (y) cos 

C A 


.(i) 


/ \ • T T X / x . II A 

w 0 g(y) sm + w 0c (y) sin 7— + 


JTX 

r 


( i ) / \ . irx (i), . irx 

w (y) sm 7- + w' (y) cos r- 

3 AC A 


(2.89) 
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Equations (2.89) are substituted into equations (2.6) and (2.7) to 
obtain the following expressions for the strains and curvatures: 


± * fl u (» 003 . fl u U) 3ln 2H „ 

x a A s A A c A 


1 ,TT\ 2 r (i) 2 (i) 2 , (i) 2 ( i) 2 , 

r(r) [w + w + (w - w ) cos . 
4 a s c sc A 


2irx 


_ (l) (l) . 2 ttx -, ^ 1 , it, 2 r (l) (l) 

- 2w w sin -—• ] + -(-) [w w„ + w w„ 
sc A 2 A s ° s c °c 


, ( l) ( l) v 2irx , ( l) 

+ (w w„ - w w. ) cos - (w w. + 
s °s c °c A s °s 


( i) , . 2 ttx -i 

w c "°c> 81n ~ ] 


(i) (i)' (i)' . 2 -itx (i)' 2 ttx 

e = v: + v sin — r- + v cos -r— + 


J{(w^ i) ') 2 + (w^’) 2 - [(w^ i} ') 2 - 


, (i)» 2, 2itx . (i)' (i)’ . 2ttx, _ 

(w ) ] cos -r— + 2w w sm } + 
c A s c A 


1r..U) 


,( i) 


(i) 


1 r \ /vo./ 

~[w w + w w„ - (w w„ - 

2 s o Q r> o r> C Oc; 


’ s c 


(i)' ’ 2ttx , (i)' ’ 

w w„ ) cos — + (w w„ + 
c °c A s °c 


(i)' ' . . 2irx, 

w w„ ) sin -r— ] 
c °s A 



81 


Y (i) = u 0 (i) ’+ u (i) 'sin^ 
xy 0 s A 


( i ) * 2 ttx 

U COS —7— + 

o X 


2tt (l) 2 ttx 2tt (l) . 2 ttx 

v cos — - — v sin -r- + 
As A A c A 


2 A sc s c sc 


(i)' (i), 2 ttx , (i) ( i) * 

W W ) COS —r— + (w W 

sc A s s 


(i) (i)' . 2ttx, 1 .it, r ( i) • ( i) • 

w w ) sin — — ] + ~(~)[w w„ - w w„ 
c c A 2 A c °s s 0 

, ( i) • (i) * , 2ttx , ( i) * 

+ (w w„ + w w„ ) cos — — + (w w„ - 
c °s s °c A s °s 


„ ( i ) ’ , . 2irx-, 1 . it , r (i) » (i) ' 

+ w w„ ) sin -r— ] + -(-)[w w„ - w w„ 
c °c A 2 A s °c c °s 


,( i) 


(i). 


+ (w w + w w„ ) cos 


2ttx 


* <W<‘V - 


) sin 2H] 
c °c A 


(l) ,TT,2 r (l) . TTX (l) TTX-, 

< = (r) [w sin — + w cos — -] 

x A s A c A 


y s 


TTX 


sin — + w 
A c 


(i) 


' TTX 

c°s J-] 


(i) 


.(i) 


- / II \ r \ A / TTX \ A ,/ . II A *1 

k - 2 (-)[vr cos — - w sin r-J 

xy A s Ac A 


.(i) 


TTX - 


(2.90) 
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The stress and moment resultants are written as 


( i )_ 

X 

x° 

(y) 

+ N (i) 
xs 

(y) 

sin 

2ttx 

X 

+ 

N (i) 
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2ttx 

+ 
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yo 
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sin 

2ttx 
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N (l) 

(y) 
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(y) 

TTX 
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M (1) 

(y) 
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TTX 
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m ( i) 
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TTX 
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(y) 
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TTX 
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xys 
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2 wx 
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where 
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a”'c- - * i(i) 2 (w^ 1)2 * »< I) 2 > . i(f) 2 (w„x"* 


.(D_ «(Di 


12 


U c 

a 


4 A ' 


, , (i) - 

’s s 


W„j i) )] + A^ ) {vi i), + J[(^ i} ) 2 + (^ i)f ) 2 ] 


' G C 


* l(w’ w (1 >'. w! w (1) )l * a' 1) [u1 1) ’. 

2 °s s °c c 26 0 


2 A s c s c °sc 


Ci) ' ’ (i) ' (On 

w„ w + w„ w - w„ w ] 

°C S °G S °S C 
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N U) = A ( ^[-(|2.)u (i) - kf) 2 (w (1) w (i) + w (i) w o + 
ys 12 X c 2X sc s°c 


.(i). 


(i)r„(i)V 1 i) i) ' ..(1)'. ' 


“o “.3> ] ‘ A 22 [v s * 5<“s 


+ w 


w„ 


+ w (i), w : )] + A^ ) Cu (i) '- V (i) + 

c °s 26 s Xc 


1 / tt \ / ( i ) (i)' , (i) (i)' ( i) • (i)' 

o (r)( w rf ~ w w + w w„ ~ w w„ 


2 X 


( i ) * ( i ) ' 

+ w UJ w„ - V )] 

s °s c °c 


„(». 4 C) C SI „»>♦ 1(I) 2 („< 1)2 - „<‘> 2 > * 


yc 


12 u x 






( i ) * (i) (i)' (i)' • ( i) , 

W W + W„ W + W„ W + W„ W + 

s c °s c °c s °c s 


w’ W (i) )] 

°s c 


c- r * j ( a ) 2(u “ >2 * »‘ 1,2 > * 1 <i> 2 <“. 3 “ 1 u 

♦ “"o“o 1>)] * ("o 11 ' 2 )] 

* i(«: «: w (i >h . 

2 °s s °c c 66 0 

2 X s c s c °sc 


( i ) T 1 (i) • (i). 

/„ w + w„ w - w„ w 
°c s °c s °s c 
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l i} - A 1 ( J ) [-(^)u (i) - kr) 2 (w (i) w (i) + w (i) w 


xys 1 6 


°c 


w (i) w 0 )] + A^ ) [v (i) '+ i(w (i) 'w (i) '+ w (i) ' w ; 

G °S 26 S 2 S G s °G 


+ )] + T 1 V (l) + 

c °s 66 s X c 


1,ir w ,( i) (i) ' (i) ( i) ’ (i)’ (i) • 

X ( 7" ) ( w w - w w + w w - w w„ 

2 A S S GC S °SG °C 


r i } i c i ) * 

+ w u V - w u V )] 

s °s c °c 


N 


(i) 


,v2 


A^^ir 1 u (l) + 7-(7) 2 (w (i) - w (i) 


.,2 


xyc 16 L X 


4V 


) + 


2V lW °ss W °c W c A 26 lV c 


1 ,, (i ) * 2 , (i)* 2 , 1, '■ ( i) • • (i)' 

5 [( “s ) - (w 0 ) ] - 5 (w os w s - w„ c w c )} 


* »66 > [ u o i) '* f v i lK ?^^ w l i>w p i),+ 

boc As 2 X S C 


( i ) ' (i) (i)' ( i ) • ' (i) 

w w + w. w + w„ w + w w + 
s c °s c °c s °c s 


w’ w (i) )] 
°s c 


M 


M 



D (l) 

{}) 

2 (i) 

w - 

D ( ^ 

(i)” 
w + 

2D 1 ^ 

{}) 

w (i) ’ 

xs 

1 1 

X 

s 

12 

s 

16 

X 

c 

,(i)_ 

D i^ ) 

{}) 

2 w (i) - 

D<« 

w U) "- 

2D^ 

(7) 

w (i) ' 

vl 

l XG 

1 1 

X 

c 

12 

c 

16 

A 

S 


D i^ 

(-) 

2 (i) 

w - 



2 D^iP 

(7) 

w (i) * 

ys 

12 

A 

s 

22 

s 

26 

A 

c 

, (i) = 

1 D iP 

(-) 

2 w (i) - 



2 D< i} 

(7) 

w (i) ' 

yc 

12 

A 

c 

22 

c 

26 

A 

s 
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M (l) = D 1 “ > (I) 2 ir (1) - D^ ) w (l) ”+ 2D^ ) (~) w (i) ’ 
xys 16 X s 26 s 66 X c 

M U) = D (i)^ ) 2 w (i)_ D (i) w (i)"_ 2 d ^)(I.) w (l) ' 

xyc 1 6 X c 26 c 66 X s 


(2.92) 


Equations (2.90) to (2.92) are used in equation (2.25), and the 
resulting equation is integrated with respect to x. For arbitrary 
,..( i) ,,.( i ) .,,( i) . , r ( i ) . „ ( i ) r, r (i) ,,,( 1 ) *„(*) ^ 

9 r> 9 9 5-W f the 

o O 5 0 5 C 

principle of stationary potential energy requires that these 
differential equations are satisfied: 

for i = 1 , 2, . . . , N 

N (l o' = 0 
xy° 

N ^ 1 N (l ^'= 0 
X xc xys 

2l N U) + N (i) '= 0 
X xs xyc 


N ( o ),= 0 
y° 

2l N (i). N (i)’ = o 
X xyc ys 

+ r/ ^ ' = 

X xys yc 


0 


(2.93) 
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V 


PC POOR 

for i ° 1 

5 ( ? )2N xi )( "o' ) * “.o’ - ( r )[( C* 

Hyo )(S o’ > *“-c ) 4">3 ), S )1 - 


! ‘*-£ tis 


0 l* 

2 xy.s 

, * c -e» 1)% ( r ,2 “s' )] - 
* [ - 6 < 2 ,, , (f) V 2 ) ] _ E * (u < ,, . w < 2 ) 


) 


v O)' = _ (f) 2 M (1). (? 2 (1 )_ 1 M (1) (1) + + 

yo X xc X x° 2 xc c °c 


H )2 »xl >( '4 ,) * “•s ) + ( ? )[(,, xy"' 

¥- 4 ' r * - 

¥' 6 o 2> '* ( r )2 “o 2)] - E * ( “o”- “o 2)) 


(2.94) 


for i = 2, 3 , . . . , N 


V (i) ’= -(I) 2 M (i) - (f) 2 (N ( ^+ ~ N (i) )(w (i) + w 0 ) + 
ys X xs X x° 2 xc s °s 

5(7) 2 N (i) (w (i) + w ) - (J)[(N (1 ]+ 

2 X xs c °c X xy° 

\ N (l) )( 3 (l ' )+ ) + 5 N (l) (S (l) + u' 0 )] - 

2 xyc c °c 2 xys s °s 
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2G*[-6^ i) '+ (f)V i} ] ~ 

S AS 




^ S AS S AS 

- E*(2w (l) - w (l_1 v/ 1 + 1 h 
S S' s 

-^ 2 <o- I “.o’ * 

„ os , . (!)[(»»)- 

I O' 6 ^ X * i O 6 o 1)+ V )] - 

2G*[-g^ l} ' + (J) 2 w^°] - 


- E*(2w (l) ~ w (l ’ 1 w (l + 1 h 

G C C 


(2.95) 


for i = N 


V 00 ’. -<i> 2 m ( n) - (I) 2 (n ( ^. « (N) )(. (,) . K„ > . 

ys X xs X x° 2 xc s °s 

1(I) 2 n (n) ( w (N) + w ) - (J)[(N (N ^+ 

2 X xs c °c X xy° 

I OX’* * 5 “la” - 
XX’* ‘XX - 


o’c-b"*' 1 ’'* (f) 2 w ( "-'>] - eV N) - 

2 s X s s s 
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'r- -w* 5 N xr><“< N) ‘ ».„> * 

„ 0s , + (I)[«m- 

I C>^ N) * »:.> * 1 O'?’* u "= )] - 


g*c- b «>' + (f) 2 w<«] - 

I C AC 


oW’ 1 ' 0 ’* (f) 2 -'"- 1 ’] - E*(w <N >- ^X-'h 
2 c AC C C 


(2.96) 


where the previously undefined functions for the ith fiber-plate are 


V (l) = -M (l) + —■ M (l) - (N ( ^~ \ N (l) )(0 (l) + w’ ) - 
ys ys A xyc y° 2 yc s °s 

5 N U) (0 (i) + w' 0 ) + (J)[(N (i ^- 5 N (i) )(w (i) 

2 ys c °c A xy° 2 xyc c 

* “• = > - 2 ^xys^s' ^ + ”• = » 

V ^ 1 ^ M (l) - (N^ + 5 N (l) )(6 (l) + w[ ) - 

yc yc A xys y° 2 yc c °c 

1 N (i) (0 (i) + Wg ) - (2)[(N (i ^+ \ N (i) )(w (i) 

2 ys s °s A xy° 2 xyc s 

* “°s> - I 0“c l)< ' V )] 

B (1) . w (1 >’ 

s s 


3 (i) = w (i) ' 
c c 


(2.97) 
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The boundary conditions for this energy formulation are 



b 

= 0 

N (1) S u (i > 

b 

0 

X yO 0 

0 

xys s 

0 


(l) 6u (l) 

b 

= 0 

N^Sv^ 

b 

0 

xyc c 

0 

yO 0 

O’ 

■ 

: (i) 6v (i) 

b 

= 0 


b 

0 

ys s 

0 

yc c 

0 



b 

= 0 


b 

= 0 

ys s 

0 

yc c 

0 



b 

= 0 

V (i) 6w (i) 

b 

= 0 

ys s 

0 

i = 1, 

yc c 

2, .... N 

0 



(2.98) 


The model remains rectangular during loading and is simply supported 
on the unloaded edges. The boundary conditions are expressed as 


at y = 0 and at y = b 
(i) ( i) . U) (!) U) (!) (i) 


U„ = U = U = V 

0 s c s 


»“>. „u>, M (1) . 0 
y° ys yc 


v ' ' = w ' ' — w ' ’=0 
CSC 


(2.99) 
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The first-order governing equations for a model with 
anisotropic fiber-plates can be expressed as a function of twenty 
unknowns for each fiber-plate. These unknowns for the ith fiber-plate 


are u^, u (l) , u (l) , v^, v (l) , v (l) , w (l) , w (l) , B (i) , B (l) , N ( ^, 
0 s c 0 s c s c s c y° •%: 


N U) , N U) , N (l >, K (1) , N (l) , M (1) , M (1) ; V (1) , and » (1) 
ys yc xy° xys xyc ys yc ys yc 


The other 


unknowns are expressed as functions of these twenty. The twenty 
differential equations for each fiber-plate are obtained by re- 
organizing equations (2.90) and (2.93) to (2.97) and are stated below: 


for i = 1 , 2, . . . , N 


1 ^ = - ^( f )^ 1 8^W^ + w B (l) - w 

" H p p °SC °r> < 


2 1 


>c s 


( i ) ^ 

» (i) - »«>) * ^Lrv 1 )- 

°c s °s c 66 (i) xy° 

A 22 


A U) A 1 A 1 u 2 

tiL. N (i). (A (i). 26 A 12 k ^ + l ( I)2 (w (i) 

.(i) y° U 16 (i) )l a 4V U s 

22 22 


2 

+ W (l) + 2w W (i) + 2w W (l) )]} 
c °s s °c c 


u s iv - r ' , o i) - 5 ( ! )(w s 1)b s 1) - 


,( i). 


,( i). 


.(i). 


- B Vi/ w 0 + w Vi 'w D - w^'w„ ) + 
c °c s °s c °c 
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A (1)2 A (1) 

U 66 (1) ’ '"xys (i) "ya 

22 22 


(a (D- A 26 A 12 }[ ,2^ (i). 1 (I) 2 ( (i) (i) 
U 16 _(i) JU c 2 a 2 ^ s W c 

A 22 


+ w^w + w (i) w 0 )]} 
s °c c °s 


u ( i} ' = - fi -/ i} - ei l) w^ l) + w 0 6 U) 


2 X ' 


°s c 


(i) 


f 


,( i ) 


(i) 


+ w 8^ /+ w 0 v^ x ' + w w' x/ ) + 

°G S °C S °S C 

A (1)2 A (1) 

,,( i) 26 .-1 , <i), 22L n (1> - 

A 66 - Tar’ 1N xyc- 7TTT V 
22 22 


^ ”“ )2 

A 22 


+ 2w w (i) - 2w 0 w^ i} )]} 
°s s °c c 


- x i)2 * 6 ' i)2 * 2 w;„ b »>) 


’c c 


( i) £ 


a; 


(i) 


,.(i) 26 r 1 f ..(i) 26 (i) 

( A 22 ~ 7 T7 N ”° — “ N -”° 


66 


y° .(i) xy° 
A 66 


A (i) A (i) u 22 

r * 

A 66 


+ 2w w< l> + 2w 0 v/ 1 *)]} 

0 3 S °C C 
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^ l} ’= - i(^ l} 6l l} + B^ l} w’ + si i} w! ) + 


2 s m g 


c c 


A (1)2 A (1) 

,.(i) 26 r 1 fM (i) A 26 M (i) 

22 (i) ys (l) xys 

06 a 66 


A k 


+ v/ * ^ w + v/ * 3 w. ) ] } 
s °c c °s 


.(!)’_ i,«( i) 2 _ o(i) 2 


= 71(8 


(i) 


, (i) 


4 K S 


+ 2w B - 2w B ) + 

C °S 3 °C C 


A <1>2 A U) 

22 (l) yc (l) xyc 

06 A 66 


' A^ 3 A^ ^ 2 2 

(i)„ *26„ A 1 6 ..pgs (i) + i ( I)2 (i) 2 _ (i) 2 

•12 (i) ;L A s 4V 1 s c 

A 66 


+ 2w v/* 3 - 2w_ w (i) )]} 
°s s °c c 


w 


w 


(i)* = fi (i) 
s s 

(i)’_ „(i) 


8, 


„ ( i) ’ J rn (i),TTv2 (i) (i) r TT. (i) M (i)-| 

B s = A T) CD 12 ( A } W S + 2D 26 ( A )8 c M ys 3 
U 22 


Q ( i) ' 1 rn (i),i,2 tl (i) (i),iK (i) M ( i) 

S o ■ 7T) CD 12 'a 1 ”o ' 2t> 26 ( A )e 3 "yo 1 
22 
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N 


N 


N 


N 


(i) ' 
y° 

U)\ 

ys 

(i)' 

yc 

(i)V 

xy° 


= 0 


(Si, N U) 

A xyc 

-(ft) N (1> 
X xys 


» (i)2 » (1 )» A (1)2 A (i) - 2A (1, A a, A U) 

N (i)' = (A (i). 12 A 66 i 6 A 22 ^ A 12 A 16 fl 26 > * 

xys U 1 1 / . \ \ , . \2 ' 

A 1 A 1 - A 1 
22 66 26 


(ft)[ft U<». i(f)V 1)2 - w <i)2 t 2w w (1) 

A A s 4 A s c °s s 


A (i) A (i) ~ A (i) A (i) 

- 2w 0 w (l) )] + ( — )[(■ — 1£_26 _ )n (i). 

° C A (i) A U) - A (i) ^ 

22 66 26 


A U) A (i) - A (i) A (i) 

/ i 6 A 22 A 12 rt 26 > (i) 1 

A (i) A (i)_A (i)2 XyC 

a 22 66 A 26 


A (i)2 A U) + A U)2 A (i) - 2A U) A (i) A (i) 

.( i) ' ,.(i) ”12 A 66 i 6 A 22 ^12 A 16 26 , , 

J = -(A,. ~ — — — ~ ■ ■■■ ) * 

XYC A U) A (i) - A U) 

A 22 A 66 A 26 


(2ZL)[-(£!)u (i) - i( 7 ) 2 (w (i) w (i) + w (i) w 0 + 
A A c 2 A sc s °c 
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A (i) A U) - A (i) A (1) 

”c "»3 >] - ^ 

A (I) A U> - A (1) 
A 22 A 66 A 26 


A 1 A 1 -- A ( 1 A 1 
,16 A 22 A 12 A 26 _UK 

( “"“T” 7^ )N xys ] 

A (l) A (l) - A 1 
A 22 66 A 26 


D (i) D (i) 

M (l) '= ~V (l) + (f^)[(2) 2 (D ( ^ } - 12 .... ,^- )w (l) - 

ys ys X X 16 (i) c 

22 


D (1) D (1) 

u 26 ..U) u 26 „(!)-, 

2< r )(D 66 jJT n s 'Jil M yc 5 " 

22 22 


(N ( ^- 1 N (l) )(g (l) + w' ) ~ o N v±; (e vx ''+w 0 ) 
y° 2 yc s °s ? n °<^ 


1 M (i)/ Q (i) ’ 

2 N ys (0 c ‘ ° c ’ 


+ (7)[(N (1 ]- \ N (i) )(w (i) + w 0 ) - 
X xy° 2 xyc c °c 

1 N (i) (w (i) + w )] 

2 xys s °s 

D (i) D U) 

m ( i), = -v (l) + (f*)[(J) 2 (D ( J } - ■ " 1 ‘ " 2 ' 7 ' - T * " ) w ( 1 } + 
yc yc X X 16 (i) s 

22 


D (i) d 1 

,,irw n (i) u 26 . fl (i) + u 26 „(!)-, 

2( r KD 66 ’ TTT )6 c + HI) M ys ] " 
22 22 


(N^^+ - N^)(g^+ w' ) - - N^(g^ + w' ) 
v "y° 2 yc Mi3 c °c J 2 ys V|3 s °s ; 

- (?)[ (N ( + \ N (i) )(w (i) + w 0 ) - 
X xy° 2 xyc s °s 


1 N (i) (w (i) + w )] 

2 xys c °c 


( 2 . 100 ) 
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for i = 1 


D 


(1 ) 


d5 1 } d5 1 } 


,,(1)' ,ir x 2 r 12 „(1 ) . „,ir w „(1 ) 12 u 26,.(1) 

V„„ = '( 7 ) L-7TT M „ Q + 2 ( rH D l6 ~ — 


ys 'X' L n d) ys 
22 


22 


J1)‘ 


jfr )( r ) M 1) ] - ( ! )2(w s 1)+ w o s ) * 


22 


S ( ’> 2 4 (1) * A (1)2 4 (,) - 2A (1 W” 

,,.( 1 ) 12 A 66 A 1 6 A 22 ^ A 12 A 1 6 A 26 x . 

UA ii p ^ 

A (1) A 0) - A (1) 

A 22 A 66 A 26 


c-r ♦ <?>‘4 , M<r ) 2 <3"i n2 * -‘ ,)2 * 


.0 ). 


( 1 ). 


6w w + 2w w )] + ( 

S °S C °G 


A (1) A (1) ~ A (1) A (1) 
A 12 A 66 A 16 A 26 

A (1) A 0) - A (1)2 
A 22 A 66 26 


) * 


A (1) A (1) - A 0) A (1) 

(N ( !> + i N (1} ) ♦ (il^ + 

7 y ° A 0) A (1) - A (1) ^ 

A 22 A 66 A 26 


~ N (1 } 
2 xyc 


* |(f) 2 (w‘' ) * K„ c ) * 


,n> 2 A<i), A (i) 2 A u)_ 2A u) A (i) A (i) 
,,„( 1 ) A 12 A 66 A 16 A 22 ^ A 12 A 16 ft 26 x , 
i' A ii ” 2 

A 0) A (1) - A (1) 

A 22 A 66 A 26 
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ORIGINAL PAGE 13 
OF POOR QUALITY 


IXfW 1 )- kf) 2 Cw^ 1 ) W (1) + W ( 1 ) W # ♦ 

X e 2 X sc s°c 


A (1) A (1) - A (1) A (1) 

w (1 )w . + ( !l„2 . „ A i6_ A l6 A 26 )n (1) + 

A (1) A (1) - X (1) ^ 

*22 *66 *26 


A (1) A (1) - A (1) A (1) 

g ' 2 - ^ 1" )0 - (?)[(«">* 

.(1).(1) .(I) 2 XyS X Xy 

*22 *66 ' *26 


1 N (1) )(B (1) 

2 xyc MP c 


N VI ')(X' + w 0 ) + i N (1) (6 (1) + w 0 )] - 
0r > 2 xys S 0 « 


„(') D «> D«) 

o*[js 7 -- 2 4 f T B (” * 

1 n ( 1 ; n (1)c n (1 ) X s 

22 22 22 


* M (2) 

- G [-22- 

2 D (1) 

D 2 

- 2 -2L b (2) + (1 

2 (2) 3 c (1 

n (2) 

°12 

d (2) 

U 22 

°22 

°22 

* Ml 
- E (w Uj - 

Q 

w (2) ) 



A S 


D ( 1) D (1) D (1) 

, r ( 1 )' , ir . 2r 12 ..( 1 ) „, irw «( 1 ) M2 U 26 ,.( 1 ) 

V yc = ’V V ~ 2 X D 1 6 " ”177" )B s 
22 22 


D 


( 1 )‘ 


v >* 


11 


22 


(D 2 (D t A o) 2 # (t). 2 A (i) A d) to 

,,.( 1 ) 12 *66 *1 6 *22 **12 *16 *26 , - 

^1 1 "* 2 

A (1) A (1) - A (1) 

22 66 26 
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A (1) A (1) ~ A (1) A (1) 

2 „ (1) „„ * 6w (1) w„ ,] . > *?? ) , 

5 3 ° ° A (1) A (1) - A 0) 

A 22 A 66 A 26 


(N „ 

y° 


(1) - 1 N (1) ) + ( 

irO O tr/* 1 V 


2 n 

yc 


i 6 A 22 A 12 A 26 (1 ) 

^-)(N - 


A (1) A (1) - A (1) ‘ 
a 22 a 66 26 




w n ) * 
°s 


{(A 


A ( 1 )2 a ( 1 ) + a (1) 2 a 0). 2A (1) A (1) A (1) 
( 1 )_ 12 A 66 A 1 6 A 22 ^ A 12 A 16 ft 26 . # 
1-4 ~ * ' 


11 


A (1) A (1) - A (1) ' 
A 22 A 66 A 26 


w + 


AC 2 X SC S-u 

a < 1) a <!>_ a (1) a <1) 

„<'>«„ )] ♦ li-JiW 1 ’.. 

C 3 A ( 1 ) A (1) A (1) ^ 

22 66 ' 26 


A (1) A (1) - A (1) A (1) 

( 16 2 ?.. IS...i£, )N (1> } ♦ (I)[(N (1) - 

A (1).(1) A (1) 2 XyS X Xy ° 

A 22 A 66 - A 26 

1 N^)(8^ + w' ) + - N^(B^+ w' )] - 

2 xyc Mi5 s »s ; 2 xys^ p c °c' J 

M (1) n (D n (1) 

G *r is- ♦ zJL. ♦ (i - Ji.)(E)V 1 )] 
G 1 L n (1) 2 n (1) e s U n (1) ; V c J 

22 22 


22 


„(2) n (2) n (2) 

G*[-i£_ + 2-21. b ( 2) + d , -li-)(S) 2 w (2) ] 

G 2 L _(1 ) 2 .(2) 6 s U n (2) } V W c J 


22 

- E*(w (1) - w (2) ) 
c c 


( 2.101 ) 
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for i = 2, 3, .... N-1 


D (1) D ( 1 , D (1) 

yO)'. _ ( S) 2 r_L2_ M U), 2 (I)( I >< 1 ). ^ 26 , Cl) 

\r o 1 *- f ■? ^ 1 «•» 1 1 6 ' 4 ' ^ 


ys L (i) ys 

22 


D 


(i) 

22 


D 


(i) 2 


* < D “ ) - nfr>^ 2 “a 1) ] - “»s> * 

u 22 


4 (l) 2 s (i) t A (i) 2 A (i)_ 2A u) A (i) A u> 
,,.(i) 12 A 66 A 1 6 A 22 ^ A 12 A 16 A 26 , * 

{(A ) * 

A (i) A (i) - A (i) 

A 22 A 66 A 26 


[ - r * ( ! >u s 1) * H )2<3 v 4 1)2 * “c 1)2 * 


.(I) 


(i). 


6w v "'w 0 + 2w' A 'w 0 )] + ( 

q Oo o 0o / - i N 


A (i) A (i) - A U) A (i) 
A 12 A 66 A 16 A 26 

A ( i ) A ( i ) — A (i)2 
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Chapter 3 


RESULTS AND DISCUSSION 


This chapter presents results from the analytical model 

described in Chapter 2 Analysis for compression-loaded laminates. 

Short-wavelength buckling results obtained using the linear analysis 

are discussed. These results are compared to the results obtained by 

Rosen [7] for unidirectional [0] -class laminates. Laminate end 

s 

shortening and stress data for short-wavelength buckling are presented 
for several multi-directional laminates. The short-wavelength 
buckling mode shapes for the laminates are described. Results 
obtained using the nonlinear analysis also are discussed. Each 
laminate has laminae with short-wavelength out-of-plane initial 
imperfections. The nonlinear behavior of a laminate is described 
using laminate end shortening, compressive load, and maximum w 
displacement. Shear stress and strain distributions within the 
laminate are presented. 

The last section of this chapter discusses failure prediction 
for compression-loaded laminates. A failure criterion is described 
and applied to the analytical results from the present model. The 
failure mode is predicted for several laminates. Approximate 
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equations for predicting the compressive failure of 0 -dominated 
laminates are presented. 

All results were obtained using the central scientific 
computing complex at NASA Langley Research Center. The results of the 
linear analysis were obtained using Cyber 170 series computers. The 
analysis of an eight-layer model required a maximum execution time of 
1000 seconds. This long execution time is the result of the number of 
iterations required by the solution algorithm to determine the minimum 
compressive stress for short-wavelength buckling and the corresponding 
half-wavelength of the mode shape. The analysis of the eight-layer 
model required approximately 39,500 decimal memory locations. The 
results of the nonlinear analysis were obtained using a Cyber 205 
vector processing computer. This analysis had a maximum of forty-one 
finite difference stations. The ratio of the maximum distance between 
adjacent stations to the plate width was 0.05 when using all forty-one 
stations. An error tolerance was specified for each unknown in the 
nonlinear analysis since the results of this analysis are approximate. 
This specified tolerance was approximately one percent of the maximum 
calculated value for each unknown. The execution time for a nonlinear 
analysis for a single load ease of a four-layer model with orthotropic 
fiber-plates was approximately twenty-five seconds. The execution 
time for a single load case of a eight-layer model with anisotropic 
fiber-plates was approximately 390 seconds. Ten incremental load 
cases were used for a typical nonlinear analysis. The analysis of 
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these four-layer and eight-layer models required approximately 1 X 10^ 
memory locations and 4 X 10^ memory locations, respectively. 

3.1 Short-Wavelength Buckling 

3.1.1 Comparison with Rosen's Results 


Rosen’s studies [7] included the short-wavelength buckling of 
the [0] g -class laminates used in structural applications. These 

laminates had fiber volume fractions greater than 0.30 and buckled 
into the shear mode shape shown in figure 1.1. He assumed that the 
laminates failed when short-wavelength buckling occurred and suggested 
that the compressive strength for these laminates was proportional to 

the matrix shear modulus as indicated in equation (1.3). The fiber's 

> 

elastic properties had a negligible effect on the predicted strength. 

The analysis used in this study can be simplified to obtain 
Rosen's results. All fiber-plates have the same w displacement for 
the shear mode shape; equation (2.55) becomes 

w^(x,y) = w s sin sin 


(3.1 ) 
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where w^ is a constant. The characteristic equation for determining 
the short-wavelength buckling of an N-lamina [0] s -class laminate is 
obtained from equations (2.56) and is expressed as 


[°n )( ?> V 2 c 66 >)( ?> 2 <5> 2+ 


4 C), A 0) 1(I) 2 t V , j2r / J ) 2 + 

( a \T 1 2 HI 4t U f V L V 


22 
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Equation (3*2) is satisfied if 


o 1 ( 1 i) (f) 4+ 2(d 1 ( ^ ) + 2d L i} )(r) 2 (S) 2+ 

„T 
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(3.3) 
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or if 
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(3.4) 


Equation (3*3) determines the critical end shortening for short- 
wavelength buckling of the ith fiber-plate when i = 2, 3. • N-1 , 
i.e., interior fiber-plates. Rosen's analysis considers interior 
fibers. The bending stiffness terms in equation (3.3) are negligible 
compared to the loading and matrix-foundation terms, and the critical 
normalized end shortening is expressed as 



2t )2 d + (£) 2 ]} 
m b 


(3.5) 


where this negative normalized end shortening corresponds to a 
compressive strain. The half-wavelength of the buckling mode shape X 

X 2 

is much smaller than the laminate width b so that the (r) term in 
equation (3.5) is negligible. The half-wavelength is discussed in 
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detail in section 3.1.1 Laminate Mode Shapes. Equations (2.1) and 
(2.2) for thicknesses and fiber volume fraction are substituted into 
equation (3.5). The constitutive equation for the ith fiber-plate 
also is substituted into equation (3.5), and the compressive stress in 
the ith fiber-plate is 


o 


(i) 

x 


V f G m 

1 - v. 


(3.6) 


The compressive stress in the composite laminate is 


a 


c 



(3.7) 


Equation (3-7) is the same as the equation (1.3) obtained by Rosen. 


The analysis used in this study both simplifies to Rosen's 
classic results and provides new results for the short-wavelength 
buckling of compression-loaded laminates. Equation (3. 4) determines 
the critical end shortening for short-wavelength buckling of the ith 
fiber-plate when i = 1 or when i = N, i.e., the outer-most fiber- 
plates in the model. Rosen's analysis does not specifically consider 
the outer-most fibers. The compressive stress in the composite 
laminate when the outer laminae buckle in a short-wavelength mode is 
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derived using the same procedure as used above. This compressive 
stress is expressed as 


o 


c 



(3.8) 


Equations (3*7) and (3.8) show the dominant role of the matrix in 
determining short-wavelength buckling response. The fiber-plates are 
so thin that their bending stiffnesses are negligible compared to the 
foundation extensional and shearing stiffnesses of the matrix. The 
compressive stress for short-wavelength buckling is a function of how 
the fiber-plate is supported. An interior fiber-plate is supported on 
two sides by matrix-foundation. The compressive stress determined by 
equation (3.8) is half the compressive stress determined by equation 
(3.7) because the outer-most fiber-plates are supported on only one 
side by matrix-foundation. 

The short-wavelength buckling of a [0] -class laminate may be 

s 

a sequence of buckling of the outer laminae. Outer-lamina buckling 
can occur when the outer-most laminae are not supported by test 
fixtures or specimen end tabs. The amplitude of the short-wavelength 
buckling mode shape may be large enough to cause matrix failure that 
leads to delamination. When one pair of outer laminae buckles and 
delaminates, the adjacent laminae become the "new" outer laminae and 
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subsequently buckle and delaminate. The failure of a [0] -class 

s 

laminate is characterized by this progressive short-wavelength 
buckling and delamination of the outer laminae. The characteristic 
"brooming" failure mode for these laminates is evidence of progressive 
buckling and delamination. 

3.1.2 Laminate Compressive Stresses and Strains 

The discussion in the previous section describes values that 

bound the compressive stress from the present model for the short- 

wavelength buckling of [0] -class laminates. The bounds for the 

s 

compressive stress result from allowing the interior laminae and the 
outer-most laminae to buckle independently. The general form for the 
present analysis requires all laminae to buckle simultaneously into 
the laminate's short-wavelength mode. The normalized compressive 
stress from the present analysis as a function of the number of 
laminae in the laminate is shown figure 3.1 . The compressive stress 
is normalized by the stress obtained from Rosen's analysis given by 
equation (3.7), and the laminate fiber volume fraction is 0.55. A 
two-lamina model focuses on outer-laminae buckling and gives the same 
value as equation (3.8) for the compressive stress from the present 
analysis. As the number of laminae in the laminate increases, the 
normalized compressive stress for short-wavelength buckling 
asymptotically approaches unity. A twenty-lamina model is similar to 
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Rosen's model and predicts a compressive stress at short-wavelength 
buckling within five percent of the compressive stress from equation 
(3.7). The significance of the behavior in figure 3.1 is that 
consistent results using the present analysis for determining the 
stress at short-wavelength buckling in different laminates require the 
laminates to have the same number of laminae. 

The laminate compressive stresses and strains for short- 
wavelength buckling of several two-lamina, four-lamina, and eight- 
lamina laminates are presented in Tables 3.1, 3-2, and 3.3, 
respectively. The critical normalized end shortening u 0 /a is equal 

to the laminate compressive strain for short-wavelength buckling. The 
half-wavelengths for the mode shapes are also included in the tables 
and are discussed in the next section on laminate mode shapes. The 
results in Table 3-1 show that these stresses are affected by laminate 
anisotropy. A combination of bending and inplane stiffnesses (see 

and h^\ respectively, in equations ( 2 . 66 )) cause the 
mn mn K J M 

compressive stress for a [+ 30 ] laminate to be thirteen percent less 

s 

than the compressive stress for a [ 0] laminate. Similar laminates 

with more than two laminae have the same differences in compressive 
stresses observed in Table 3.1. The compressive stress for a [90] 

O 

laminate is approximately the same as the compressive stress for a 

[0] laminate. These two specially-orthotropic laminates demonstrate 
s 
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that matrix-foundation components dominate the laminate compressive 
stress. All compressive stresses for the four-lamina laminates (Table 

T 

3.2) are approximately the same. These laminates do not have the A r 

1 o 

T 

and A 2 g inplane stiffness terms that affect the compressive stress 

of some two-lamina laminates. The matrix-foundation components 
dominate the compressive stresses for the four-lamina laminates. The 
compressive stresses in Table 3.3 for the [0^] , C ( ±-45 ) ] 3 » C + 45/0/— 

45/90] , and [90/+45/0/-45] laminates are approximately the same, 
s s 

The matrix-foundation components dominate the compressive stress for 

these laminates. The compressive stresses for the [0/±45/90] and 

s 

[0/+45/90/-45] laminates are thirteen percent lower than the 
s 

compressive stress for the [0 i| ] s laminate. These lower stresses may 

0 

be due to the laminates having 0 laminae as the outer-most laminae 

O 

and not having any interior 0 laminae. 

The effect of fiber volume fraction on the normalized 
compressive stress for short-wavelength buckling is shown in figure 
3.2. The compressive stress is normalized by the compressive stress 
for short-wavelength buckling of a laminate with fiber volume fraction 
equal to 0.55. The normalized compressive stress is determined for 
0.25 S Vj, £ 0.85 which includes most all composite laminates in 

structural applications. These results agree with similar results by 
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Rosen [7]. The results from the present analysis show that the curve 
in figure 3.2 determines the compressive stress for short-wavelength 
buckling of a laminate with any stacking sequence and any thickness. 

3.1.3 Laminate Mode Shapes 

The laminate mode shape at short-wavelength buckling has half- 
waves that are either normal or skewed to the direction of loading. 
Also, the shear mode (see figure 1.1) is predicted for all the 
symmetric laminates in this study. The mode shapes for several 
laminates are shown in figures 3.3 to 3.9. The skew angle is denoted 

as The mode shape for the [0] laminate has normal waves and is 

s 

illustrated in figure 3.3* A similar mode shape is predicted for all 

balanced, symmetric laminates. The only difference among these mode 

shapes is the half-wavelength A. The mode shapes for the [+10] , 

s 

[+20] , [+30] , [+45] , [+60] , and [+80] laminates have skewed waves 
s s s s s 

and are illustrated in figures 3-4 to 3.9, respectively. The skew 

angle is plotted as a function of the fiber orientation in figure 3.10 

for these laminates. Generally, the skew angle is not the same as the 

fiber orientation. The skew angle increases to a maximum at 

O 

approximately 0 = 45 . The laminate mode shape at short-wavelength 
buckling is affected by the fiber-plate bending stiffnesses in 
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contrast to the compressive stress which is unaffected by these 
bending stiffnesses. The bending stiffness contribution to the mode 
shape is discussed in the next paragraphs. 

The half-wavelength of the mode shape for a [0] -class 

3 

laminate is determined exactly using equations (3.3) and (3.4). The 
compressive strain (normalized end shortening) in these equations is 
minimized with respect to the half-wavelength X to obtain 


X 



-1 


[D 22 ^ 


°1 ^m ,tt. 2-, i -1 /H 

2(1-V f rb ; JJ 


(3.9) 


where c 


(i) 

1 


is defined as (see equation (2.6)) 


ej 1} = 1, outer-lamina buckling 

= 2, inner-lamina buckling (3.10) 

The half-wavelength of the mode shape for any laminate can be 
determined approximately by numerically minimizing the compressive 
strain with respect to X. The compressive strain as a function of X 
is shown in figure 3.11 for several laminates. The compressive strain 
for the four-lamina laminates decreases rapidly to a minimum and 
increases very gradually as the half-wavelength increases. The 
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minimum compressive strains are observed easily by expanding the 
ordinate in the region of interest and re-plotting the results as 
shown in figure 3.12 for a [0 o ] laminate. The compressive strain for 

the C 0/1-45/90] s laminate in figure 3.11 decreases rapidly to a minimum 

and then increases rapidly in the neighborhood of the minimum as the 
half-wavelength increases. The minimum compressive strain for this 
laminate is clearly shown in figure 3.11. The half-wavelengths for 
two-, four-, and eight-lamina laminates are presented in Tables 3.1, 
3.2, and 3.3* respectively. 

The half-wavelength as a function of fiber orientation 0 is 

shown in figures 3.13 and 3.14 for [+e] laminates and [±0] 

s s 

o o 

laminates, respectively. The results in figure 3.13 for 0 S 0 £ 90 
show that the half-wavelength increases slightly and decreases. The 
slight increase may be the effect of anisotropy on laminate behavior 
discussed previously for compressive stress. The decrease is similar 
to the behavior of 1 as a function of 9. Equation (3.9) shows 

that D is used to calculate the half-wavelength. The results in 

O O 

figure 3.14 for 0 < 0 S 90 show that the half-wavelength 

decreases. The [ ±0 ] laminate behavior is not affected by the inplane 

s 

anisotropic terms that influence the [+0] laminate behavior. The 

s 

behavior in figure 3.14 may also reflect the influence of 


on the 
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half-wavelength. The half-wavelength as a function of the laminate 
fiber volume fraction is shown in figure 3.15. The response for the 
[ 0 ] and [±45] laminates is used in the figure as typical behavior. 

c. S S 

The half-wavelength increases with increasing fiber volume fraction 
for both laminates. Specifically, the half-wavelength for these 
laminates with V^, = 0.85 is approximately twice the half-wavelength 

for the corresponding laminates with V = 0.25. 

The half-wavelength results of the present study may explain 
the scatter in compressive strength data for composite laminates. The 
range of half-wavelengths for the buckling of laminates with 
V = 0.55 in this study is from 0.0140 inches to 0.2255 inches. Most 

of the half -wavelengths are on the order of 0(10 1 ) inches. Some 
compression test methods for composite laminates use specimens with 
short test sections (e.g., 0.50 in., [61]) or require elaborate 
fixtures that completely support the test specimen along the length 
[62]. These methods may inhibit the natural failure mode for a 
laminate that deforms in a short-wavelength buckling mode shape. 

These methods may cause the half-wavelength to be less than the half- 
wavelength corresponding to the minimum compressive strain if short- 
wavelength buckling ever occurs. The results in figure 3.11 show that 
significantly higher compressive strains can occur if the half- 
wavelength is 0.1 inch less than the half-wavelength corresponding to 
the minimum compressive strain. 
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3.2 Nonlinear Behavior 

The nonlinear behavior of four laminates is presented in this 

section. The laminate nonlinear behavior is the result of short- 

wavelength out-of-plane initial imperfections in each lamina. The 

laminate orientations are [0_] . [0/90] , [+115] , and [+45/0/-45/90] . 

d. s s s s 

These laminates illustrate the range of application for the present 

analysis. Specifically, these laminates illustrate: 


[ °2 ] s 

Orthotropic 

laminae, 

unidirectional laminate 

[0/90] s 

Orthotropic 

laminae, 

multi -directional 


laminate 



[ ± 45] s 

Anisotropic 

laminae, 

multi -directional 


laminate 



[+45/0/-45/90] 

s 

Orthotropic 

and anisotropic laminae, multi- 


directional 

laminate 



The nonlinear laminate behavior is discussed for two imperfection- 
amplitude-to-lamina-thickness ratios w 0 /t. These ratios are 

w 0 /t = 0.1 and w 0 /t = 0.5. The imperfection amplitudes for 

anisotropic fiber-plates, w 0 and w 0 , are simplified by letting 

s c 
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w 0 = w 0 = w 0 . All fiber-plates have the same imperfection. The shape 
s c 

of the imperfection is the same as the short-wavelength buckling mode 
shape for the laminate. 

3.2.1 Displacements and Inplane Stresses 

The nonlinear behavior of a [0_3 laminate with initial 

2 s 

imperfections is shown in figure 3.16. The ordinate is the laminate 
compressive stress normalized by the laminate compressive stress for 
short-wavelength buckling. The abscissa is the laminate end 
shortening normalized by the laminate end shortening for short- 
wavelength buckling. The laminate response for w 0 /t = 0. (no initial 

imperfections) is included for comparison with the laminate response 
for w 0 /t = 0.1 and 0.5. The results in the figure show that this 

laminate behaves like a wide column, i.e., the slope of each response 
curve (laminate stiffness) approaches zero as the end shortening is 
increased. This change in stiffness is dramatic for the laminate with 

w 0 /t = 0. but is more gradual for the laminates with w 0 /t = 0.1 and 

0.5. A maximum compressive stress that corresponds to the maximum 
load-carrying capacity of the laminate is defined by this wide-column 
response. The maximum stresses for laminates with w 0 /t = 0.1 and 0.5 
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are approximately five and twenty-five percent lower, respectively, 
than the maximum stress for the laminate with w 0 /t = 0. 

O 

The stress resultant for a 0 lamina is plotted in 

figure 3.17 for u 0 /(u 0 ) cr = 0.5, 0.8, and.1,0. The laminate end 
shortening for short-wavelength buckling is (u 0 ) . The initial 

imperfection for the [0 o ] laminate is w 0 /t = 0.5. The results in 

the figure are typical of a plate with an out-of-plane initial 
imperfection. The stress resultant at the edges of the lamina 
(y/b = 0. and y/b = 1.) is greater than the stress resultant at the 
center of the lamina (y/b = 0.5). The difference between the stress 
resultants for these locations increases as the laminate end 
shortening increases. The stress resultant at the edges is more than 
twenty-five percent greater than the stress resultant at the center 
when u 0 /(u 0 ) Qr = 1.0. 

O 

The w displacement for the outer-most 0 lamina along 
x = X/2 is plotted in figure 3.18 for u 0 /(u 0 ) cr = 0.5, 0.8, and 1.0. 

The initial imperfection for this laminate is w 0 /t = 0.5. Again, the 

results in the figure are typical for a plate with an out-of-plane 
initial imperfection. The shape of this w-displacement curve is 
approximately parabolic for u 0 /(u 0 ) cr ,= 0.5. This shape changes as 
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u 0 /( u 0 ).„ increases. For u 0 /(u 0 ) = 1.0, the magnitude of the 

slope, w »y^» i- 3 very large near the lamina edges and the w 
displacement is approximately constant for 0.3 £ y/b £ 0.7. 

The maximum w displacement as a 'function of end shortening 
for a [0-] laminate is shown in figure 3.19. The laminate response 

C b 

for w 0 /t = 0. is included for comparison with the laminate response 

for w 0 /t = 0.1 and 0.5. The maximum w displacement occurs at 

y = b/2 as illustrated in figure 3.18. The results in figure 3.19 
show that significantly large w displacements occur for the 
imperfect laminates. The w displacement is approximately equal to a 
lamina thickness for the laminate with w 0 /t = 0.1 for 

u 0 /(u 0 ) cr = 1.0 and almost equal to two lamina thicknesses for the 

laminate with w o /t = 0.5 for u 0 /(u 0 ) cr = 1.0. These large w 

displacements produce large interlaminar shear stresses as discussed 
in section 3.2.2 Interlaminar Strains. 

The nonlinear behavior for a [0/90] laminate with initial 

s 

imperfections is shown in figures 3.20 and 3.21. Laminate compressive 

stress as a function of laminate end shortening is plotted in figure 

3.20. This multi-directional laminate has the same wide-column 

behavior as the unidirectional laminate (cf., figure 3.16). The 

maximum stresses for [0/90] laminates with w„/t = 0.1 and 0.5 are 

s ° 
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approximately seven and twenty-eight percent lower, respectively, than 
the maximum stress for this laminate with w 0 /t = 0. These 

differences are similar to those observed for the corresponding [o _] 

£_ S 

laminates. The maximum w displacement as a function of end 

shortening for a [0/90] laminate is shown in figure 3.21. The 

s 

results in this figure for the multi-directional laminate are 
approximately the same as the results in figure 3.19 for the 
unidirectional laminate. Generally, the results for the [0/90] 

s 

laminates are approximately the same as the corresponding results for 
the [0J laminates. 

c. S 

The nonlinear behavior for a [±45] laminate with initial 

s 

imperfections is shown in figure 3.22. The normalized laminate 
compressive stress is plotted as a function of normalized laminate end 
shortening for u 0 /(u 0 ) cr S 1.0. Results for u 0 /(u 0 ) cr ,> 1.0 did not 

converge to within the specified tolerances. The mode shape for 
short-wavelength buckling of this laminate changes from having normal 
waves to having skewed waves for u 0 /(u 0 ) cr > 1.0. The lack of 

convergence is caused by the changing mode shape. The results in the 
figure approach the wide-column behavior discussed previously. The 
maximum compressive stresses in the figure for [±45] laminates with 

w 0 /t = 0.1 and 0.5 are approximately five and twelve percent lower, 

respectively, than the compressive stress for short-wavelength 
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buckling of this laminate with w 0 /t = 0. The results for the [±45] 

3 

laminate with w 0 /t = 0.1 are approximately the same as the results 

for similar [ 0 „] and [0/90] laminates. The results for the [±45] 
c. s s s 

laminate with w 0 /t = 0.5 suggest that this imperfection has less 
effect on the compressive stress of a [±45] laminate than on the 

3 

compressive stresses of [0 2 ] g and C 0/90 ] s laminates. 

O 

The N and N stress resultants for a 45 lamina are 
x xy 

plotted in figure 3.23 for u 0 /(u 0 ) c ^= 0.25 and 0.50. The initial 

imperfection for this [±45] laminate is w /t = 0.5. Laminates with 

s ° 

u o /( u o) cr > 0.50 have compressive stresses more than an order of 
magnitude greater than the typical compressive strength for a [±45] 

s 

laminate, and results for these laminates are not included in the 
figure. The maximum N occurs at the laminate edges although N 

X X 

is approximately constant across the laminate width. The stress 
resultant N is constant across the laminate width. The stress 

xy 

resultant N is more than twice N for both values of u„/(u„) 

x xy 0 0 cr 

however, N is sufficient to initiate failure within the laminate, 
xy 

and this failure is discussed in the section on laminate failure 


predictions 
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Results for the w displacements in the [±45] laminates are 

s 

presented in figures 3.24 and 3.25. The w displacement for the 

O 

outer-most 45 lamina along x * X/2 is plotted in figure 3.24 for 
u 0 /( u 0 ) cr = 0*25 and 0.50. The initial imperfection for this laminate 

is w 0 /t = 0.5. The shape of these w-displacement curves is 

approximately parabolic for both values of u 0 /(u 0 ) cr> . The maximum w 

displacement as a function of end shortening is shown in figure 3.25. 
The laminate response for w 0 /t = 0. is coincident with the ordinate. 

The laminate response for w 0 /t = 0.1 and for w 0 /t = 0.5 are 

presented also. Significantly large w displacements occur for the 
imperfect laminates. The maximum w displacement is greater than 1.5 
lamina thicknesses for the laminate with w 0 /t = 0.1 for 

u 0 /( u o) cp = 1*0 and is approximately equal to two lamina thicknesses 

for the laminate with w„/t = 0.5 for u„/(u.) = 0.9. These w 

o o o C ^ 

displacements are larger than the corresponding w displacements for 
the [0„] and [0/90] laminates because the axial bending stiffness 

<L S S 

for the [±45] laminates is less than the axial bending stiffness for 
s 

o 

the 0 -dominated laminates. 

The nonlinear behavior for a [+45/0/-45/90] laminate with 

s 

initial imperfections is shown in figure 3.26. The normalized 
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laminate compressive stress is plotted as a function of normalized 
laminate end shortening for u 0 /(u 0 ) cr> £ 1.0. Results for 

u 0 / ( u 0 ) Qp > 1*0 did not converge to within the specified tolerances. 

The lack of convergence is caused by changes in the laminate short- 

wavelength mode shape and is similar to the lack of convergence noted 

for the [±45l results. The results in the figure approach the wide- 
s 

column behavior discussed previously. The maximum compressive 

stresses in the figure for [+45/0/-45/90] laminates with w 0 /t * 0.1 

s 

and 0.5 are approximately seventeen and thirty-six percent lower, 
respectively than the compressive stress for short-wavelength buckling 
of this laminate with w 0 /t = 0. These differences for an eight- 

lamina laminate are larger than the corresponding differences for the 
four-lamina laminates and are related to laminate thickness. 

O 

The N and N stress resultants for a 45 lamina are 
x xy 

plotted for u 0 /(u 0 ) c ^= 0.25 in figure 3.27 and for u 0 /(u 0 ) cr = 0.50 

O 

in figure 3.28. The corresponding N x values for a 0 lamina also 
are plotted in the figures. The initial imperfection for the laminate 

_ O O 

is w„/t = 0.5. The maximum N in the 0 and 45 laminae occurs at 
0 x 

O 

the laminate edges. For the 45 lamina N x is approximately constant 

across the laminate width, and N is constant across the laminate 

xy 
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width. The difference between N for the 0 lamina and N for the 

x x 

0 

45 lamina is approximately the same as the difference between the 
corresponding axial stiffnesses. 

Results for the w displacements in the [+45/0/-45/90] 

s 

laminates are presented in figures 3.29 and 3.30. The w 

0 _ 

displacement for the outer-most lamina (45 lamina) along x = \/2 is 
plotted in figure 3.29 for u 0 /(u 0 ) cr = 0.25 and 0.50. The initial 

imperfection for this laminate is w 0 /t = 0.5. The shape of these w- 

displacement curves is approximately parabolic for both values of 
u 0 /(u 0 ) or . The maximum w displacement as a function of end 

shortening is shown in figure 3.30. The laminate response for 

w 0 /t =0. is coincident with the ordinate. The laminate response for 

w 0 /t = 0.1 and for w Q /t = 0.5 are presented also. The results are 

plotted using the same scale as that in figures 3.19, 3.21, and 3.25 
and indicate that the w displacements for the eight-lamina laminate 
are much smaller than the w displacements for the four-lamina 
laminates. The w displacements for the eight-lamina laminate are 
only one-quarter to one-third the w displacements for the four- 
lamina laminates. The w displacements for the eight-lamina laminate 
are smaller than the w displacements for the four-lamina laminates 
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because the axial bending stiffness for the eight-lamina laminate is 
greater than the axial bending stiffness for the four-lamina laminate. 

3.2.2 Interlaminar Strains 

The strains in each matrix-foundation are expressed in 
equations (2.20) as functions of the w displacements of the adjacent 
fiber-plates. These strains are interlaminar normal and shearing 

strains. The interlaminar normal strain is negligible for all 

the cases in this study since the deformation of every laminate is 
dominated by the shear mode (see figure 1.1). This deformation also 

causes the interlaminar shear strains Y^ and Y^ to be 

myz mxz 

constant within each matrix-foundation region. These strains are 

largest at the interface between the two outer laminae because the 

gradients of the w displacements are largest at this interface. The 

interlaminar shear strains are represented by Y and Y and are 

yz xz 

reported for the interface between the two outer laminae in subsequent 

discussion. The maximum Y occurs along x = mX/2, m = 1 , 2, 

yz 

3, ... (see figure 2.5a) since the w , component to this shear 

strain has a maximum along these lines. The maximum Y occurs 

XZ 

along nodal lines, x = nX , n = 0, 1 , 2, ... (see figure 2.5a) since 
the w,^ component to this shear strain has a maximum along these 
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lines. These interlaminar strains are caused by the geometrically 
nonlinear behavior of the laminae and are different from the 
interlaminar strains that result from lamina material property 
differences [ 63 ]. 

Interlaminar shear strain distributions in a [0_] laminate 

2 s 

are shown in figures 3-31, 3-32, and 3.33. The strains in figures 
3.31 and 3.32 are for u 0 /(u 0 ) cr = 0.30 and 0.60, respectively, and for 

w 0 /t = 0.1. The results in the figures show that the magnitude of 

Y is a maximum at the laminate edges. The maximum and minimum 

yz 

values for w,^ occur at y/b = 0. and 1., respectively. The results 

in the figures also show that Y has a parabolic distribution 

xz 

across the laminate width and is a maximum at the laminate center. 

The maximum value for w,^ occurs at y/b = 0.5. The strains in figure 

3.33 are for u Q /(u 0 ) cr = 0.30 and for w 0 /t = 0.5. The behavior for 

Y and Y for w„/t = 0.5 is approximately the same as the 

yz xz 0 

behavior for Y and Y for w„/t = 0.1; however, the maximum 
yz xz 0 

Y and the maximum Y at w„/t = 0.5 are greater than the maximum 

yz xz 0 

Y and the maximum Y , respectively, for both load cases at 

yz xz ’ ^ 

w 0 /t = 0.1. The w 0 /t = 0.5 imperfection causes interlaminar strains 

for u /(uj = 0.3 that are greater than the interlaminar strains 
0 0 cr 
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for w/t - 0.1 at twice the load level. The maximum Y 

0 xz 

approaches four percent at an end shortening that is only thirty 

percent of the end shortening at short-wavelength buckling. The 

results in figures 3-31 to 3.33 also illustrate that the maximum Y 

X z 


is more than an order of magnitude greater- than the maximum 



Results for Y xz only are presented subsequently. 


0 0 

The Y distribution at the 0 /90 interface in a [0/90] 
xz s 

O O 

laminate is shown in figure 3.34 and at the +45 /-45 interface in a 

[±45] laminate is shown in figure 3-35. Both laminates are loaded to 
s 

U 0 /(u 0 ) cr = 0.20, and the results in the figures are similar. These 

results are also similar to the Y results in figures 3-31 and 3-33 

X z 

for a [0 ] laminate with u 0 /(u 0 ) = 0.30. 

C. S Ci 


0 o 

The Y distribution at the +45 /0 interface in a 
xz 

[+45/0/-45/90] laminate is shown in figure 3.36. The laminate is 
s 

loaded to u„/(u„) = 0.10. The results in the figure are similar to 

O O Qp 

those reported for the four-lamina laminates except that this laminate 
is loaded only to u 0 /(u 0 ) cr = 0.10. The w >x component to Y xz for 

the [+45/0/-45/90] laminate is large because the wavelength of the 
s 

mode shape for this laminate is extremely small, i.e., X = 0.0477 in. 
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(Table 3.3). The w-gradients for this laminate are large although the 
w displacements are small (figure 3-30) when compared to similar 
results for the four-lamina laminates. 

3.2.3 Effects of Fiber Volume Fraction 

The fiber volume fraction V f affects the short-wavelength 

buckling behavior of compression-loaded laminates. The effect of V 

on the compressive stress at short-wavelength buckling and on the 
half-wavelength of the buckling mode are illustrated in figures 3.2 
and 3.15, respectively. 

The fiber volume fraction also affects the nonlinear behavior 
of compression-loaded laminates with short-wavelength imperfections. 
The nonlinear results reported in the preceeding paragraphs are for 
laminates with V f = 0.55. This section examines the behavior of 

laminates with V f = 0.45 and with V f = 0.65. These V^, bound 

the typical laminates with structural applications. Typical results 
are presented in the following paragraphs for [0 ] and [±45l 

c S 3 

laminates. 

The effect of V f on the behavior of C ° 2 1 s laminates is shown 

in figures 3-37 to 3.39. Laminate compressive stress as a function of 
the laminate end shortening is presented in figure 3.37. The laminate 
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response for w 0 /t = 0. is included for comparison with the laminate 

response for w 0 /t = 0.1 and 0.5. The results for V^= 0.H5 and 

V f = 0.65 in this figure bound the results for V f = 0.55 in figure 

3.16. The laminates with V^= 0.65 are the stiffest, as expected. The 

effect of fiber volume fraction on laminate compressive-stress versus 
end-shortening response is more significant for w 0 /t = 0.5 than for 

w 0 /t =0.1. The compressive stress for V^= 0.65 with 

u o/(u 0 ) cr = is more than ten percent greater than the compressive 

stress for V f = 0.M5 with u 0 /(u 0 ) cr = 1.0 when w 0 /t = 0.5. Similar 

results differ by two percent when w o /t = 0.1. The maximum w 

displacement as a function of laminate end shortening is shown in 
figure 3.38. These results also bound the results for V f = 0.55 in 

figure 3.19. Significantly large w displacements are observed for 
this laminate at each fiber volume fraction. The w displacements 
for V f = 0.115 are greater than or equal to a lamina thickness for 

w 0 /t = 0.1. The w displacements for = 0.65 are greater than or 

equal to 1.5 lamina thicknesses for w 0 /t = 0.5. The laminates with 

V^= 0.65 have the largest w displacements of the laminates studied 

because the matrix-foundation regions are thinnest for V^= 0.65. 

This thin region may allow large w displacements for the fiber- 
plates. The interlaminar shear strain distribution across the 
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laminate is plotted in figure 3.39 for u 0 /(u 0 ) cr< = 0.30. The 

corresponding results for V f = 0.55 are plotted in figures 3.31 and 

3.33* The effect of fiber volume fraction on y is more 

xz 

significant for w./t = 0.5 than for w_/t =0.1. The Y for 
0 0 xz 

V = 0.65 at y/b =0.5 is more than forty-three percent greater than 

Y xz for V^= 0.45 at y/b = 0.5 when w 0 /t = 0.5. Similar results 

differ by twenty-eight percent when w 0 /t = 0.1. The large w 

displacements shown in figure 3-38 may lead to large w-gradients that 
contribute to Y 

xz 

The effect of V„ on the behavior of [+45] laminates is shown 
t s 

in figures 3.40 to 3.42. Laminate compressive stress as a function of 
the laminate end shortening is presented in figure 3.40. The results 
in this figure for V f = 0.45 and V f = 0.65 are within five percent 

of the results in figure 3.22 for V^= 0.55. The laminates with 

Vj,= 0.65 are the stiffest. The effect of fiber volume fraction on 

the results in figure 3*^0 is noticeable but not significant because a 
change in V^. has a noticeable but not significant change in the 

axial extensional stiffness for this laminate. The maximum w 
displacement as a function of laminate end shortening is shown in 
figure 3*41. The laminate response for w 0 /t = 0. is coincident with 


the ordinate. These results are also very similar to the results in 
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figure 3.25 for V = 0.55. Significantly large w displacements are 

observed for this laminate at each fiber volume fraction, and the 
largest w displacements occur in laminates with V f = 0.65. The 

interlaminar shear strain distribution across the laminate at the 

+45 /-45 interface is plotted in figure 3.42 for u o /(u 0 ) cr = 0 . 30 . 

The corresponding results for V = 0.55 are plotted in figure 3.35. 

The Y xz for V f = 0.65 at y/b = 0.5 is more than thirty percent 

greater than the Y for V = 0.65 at y/b = 0.5 when w 0 /t = 0.5. 

x z r 

Similar results are observed when w 0 /t = 0.1. The large w 

displacements shown in figure 3-41 may lead to large w-displacement 
gradients that contribute to Y 

X z 

3.3 Laminate Failure Predictions 


3.3.1 Dominant Mechanisms 

This section applies the results of the present analysis for 
the short-wave length buckling response and for the geometrically 
nonlinear response of a laminate to the failure of compression- loaded 
laminates. Specifically, laminate failure initiated by outer-lamina 
buckling, by interlaminar shear strains from lamina imperfections, or 
by inplane shearing stresses is considered. Outer-lamina buckling 
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occurs at an outer-lamina stress determined by equation (3.8). This 
phenomenon is independent of the fiber orientation of the outer 
laminae since short-wavelength buckling behavior is dominated by the 
matrix contributions. This phenomenon is not observed for many 
laminates because other mechanisms dominate their failure. An example 
of such a mechanism is interlaminar shearing caused by lamina short- 
wavelength imperfections. A simple maximum-shear-strain criterion is 
used in this study to predict failure due to this mechanism. The 
shear strain for matrix failure is used for this maximum shear strain. 
The nonlinear analysis is used to calculate the laminate compressive 
stress for this interlaminar shearing failure. Another mechanism that 
may dominate laminate failure is the inplane shearing at the fiber- 
matrix interface and in the epoxy matrix between fibers. This inplane 
shearing has been referred to as matrix shearing and has been shown to 
initiate failure in compression-loaded [±45] s -class laminates [64], A 

simple maximum-shear-stress criterion is used in this study to predict 
failure due to matrix shearing. The nonlinear analysis derived herein 
also is used to calculate the laminate compressive stress for inplane 
shear failure. 

The compressive strength a c is shown as a function of lamina 

orientation in figure 3.^3 for [ ± 0 ] laminates. The laminate 

s 

compressive strength is normalized by the compressive stress for 
short-wavelength buckling of a [0„] laminate, i.e., 0 O = 433-5 ksi 

<_ 5 
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(Table 3.2). Short-wavelength buckling of the outer laminae occurs at 
o / a 0 = 0.66 (289.6 ksi) for all 0 and is indicated by the horizontal 

line in the figure. The axial compressive stress for laminate failure 
by short-wavelength buckling of the outer laminae is less than five 
percent greater than the compressive stress for interlaminar shear 

0 O 

failure when 0 S 0 S 15 . Laminate failure due to interlaminar 
shearing for laminates with w 0 /t = 0.1 and for laminates with 

w 0 /t =0.5 is plotted in the figure for a maximum shear strain 

(Y ) = 0.036 [54]. The laminate compressive strength due to 

xz max 

interlaminar shearing decreases as 0 increases. Laminate failure 

due to matrix shearing is plotted in the figure for a maximum shear 

stress in the principal material coordinate system (t,„) = 13.8 ksi 

12 max 

[64]. The laminate compressive strength due to matrix shearing 

0 0 

approaches infinity for 0 near 0 and for 0 near 90 . Results for 

O 0 

inplane matrix shearing are plotted for 10 £ 0 5 85 . The results 

in figure 3.43 suggest that laminate failure for [±0] laminates may 

s 

be due to short-wavelength buckling of the outer laminae for 

O O — 

0 £ 0 < 15 when w 0 /t <0.1. The initial imperfections w 0 /t < 0.1 

are very small and may rarely occur in typical laminates. Laminate 

failure for [±0] laminates is due to interlaminar shearing for 
s 

0 0 

0 S 0 < 15 when w 0 /t £ 0.1. The compressive strength due to 
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interlaminar shearing is a function of w 0 /t. The results in the 
figure also suggest that laminate failure is due to inplane matrix 

O O 

shearing for 15 £ 0 S 75 . The compressive strength due to matrix 
shearing is not a function of w 0 /t. Results from the present model 

O O 

for 75 < 9 £ 90 are not applicable to laminate failure. Typical 

failure in these laminates appears to be initiated by compressive 
failure of the matrix. The present model treats the matrix as an 
elastic foundation and axial loading of the matrix-foundation is not 
considered. Nevertheless, the present model predicts laminate 
failures initiated by outer-lamina buckling, by interlaminar shearing, 

O 0 

or by inplane matrix shearing for 0 S 0 S 75 . The present model is 
unique in its ability to predict the compressive strength as a 
function of short-wavelength buckling and shear failures for such a 
variety of laminates. 


O 

3.3.2 Simple Equations for 0 -Dominated Laminates 


Many laminates used in structural applications have 0 
laminae, and the behavior of these laminae often dominates the 

O 

behavior of the laminate. Such laminates are referred to as 0 - 

0 

dominated laminates. Some examples of 0 -dominated laminates are 


[0/90] and [0/±45/90] . The compressive failure of a 0 -dominated 
s s 
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laminate can be predicted by determining the failure of the 0 

laminae. The present analysis can be used to determine the 

0 

compressive strength of 0 laminae, and the results in figure 3.43 

show that many [0] -class laminates fail due to interlaminar shearing, 
s 

o 

This failure mode can cause failure of 0 laminae in a laminate. 

A simple method can be used to predict the compressive failure 

O 

of 0 -dominated laminates. The method is referred to herein as the 
stiffness-ratio method and is outlined as follows: 


O 

1 . Determine the compressive strength of a 0 lamina by using 
the results from the present analysis in figure 3* ^3; 

0 

2. Calculate the load in the 0 laminae of the given laminate 

O 

when the 0 laminae will fail; 

3. Calculate the load in all other lamina of the given 

0 

laminate when the 0 laminae will fail by using 



(3.11 ) 


where 
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P = load in a 0-oriented lamina 

O 

0 

P 0 = load in a 0 lamina 

E = Young's modulus in the x-direction for a 6-oriented 
xe 

lamina (see figure 2.4) 

E „ is determined from [12] 

X0 



1 4 ,1 

TT- cos 0 + (- rr- 

E 11 G 12 E 11 


12 . 2 . 2 . 
— )cos 0 sin 0 


+ 7T— sin 4 © 
E 22 


( 3 . 12 ) 


where E^, E^, G^, and 2 are lamina properties in 

the principal material coordinate system; 

4. Sum the loads of all laminae to determine the laminate 
failure load; 

5. Calculate the laminate compressive strength. 

The stiffness-ratio method assumes that the laminate is loaded by 
uniform end shortening (i.e., constant strain) and that the laminate 


has linear stress-strain behavior 
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The stiffness-ratio method results in a simple equation for 

0 

the compressive strength of a 0 -dominated laminate 


a lam 



v 

L 

i=l 


,(i) 

J xe 


(3.13) 


where 


°lam = GOm P ress i ve strength of the laminate 

O 

a Q = compressive strength of a 0 lamina 
N = number of laminae in the laminate 

The predicted laminate strength using equation (3.13) is compared to 

O 

the experimental strength for three 0 -dominated laminates in Table 
3.4. All laminates have V^.= 0.55, and the compressive strength of a 

O 

0 lamina is determined using figure 3.43 with o /a = 0.45. The 

0 

experimental strengths are typical results from extensive testing. 

The agreement between predicted strength and experimental strength for 
this limited number of laminates is excellent. The maximum difference 
between the predicted strengths and the experimental strengths is less 
than five percent. 



Chapter 4 


CONCLUSIONS 


4.1 Concluding Remarks 

This investigation studies the short-wavelength buckling (or 
microbuckling) of multi-directional composite laminates loaded in 
uniaxial compression. This investigation also studies the 
interlaminar shear failures due to short-wavelength initial 
imperfections for the laminae and the inplane shear failures in these 
laminates. A laminate model is presented that idealizes each lamina. 
The fibers in the lamina are modeled as a plate, and the matrix in the 
lamina is modeled as an elastic foundation. The model is applied to 
symmetric laminates having linear material behavior. The laminates 
are loaded in uniform end shortening and are simply supported on all 
edges. 


The present model is used to determine linear and nonlinear 
laminate responses. A linear analysis is derived to determine the 
short-wavelength buckling response of composite laminates. The out- 
of-plane w displacement for each plate is expressed as a 
trigonometric series in the half-wavelength of the mode shape for 
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laminate short-wavelength buckling. Results from this linear analysis 
are also compared to previous results for unidirectional laminates. 

The present linear analysis for laminates with no initial 
imperfections is generalized to obtain a nonlinear analysis for the 
response of laminates with short-wavelength initial imperfections. 

This nonlinear analysis is derived using nonlinear strain-displacement 
relations. The results of the present linear and nonlinear analyses 
are used to develop a compressive failure criterion for composite 
laminates. 

The present linear analysis is used to determine the laminate 
stresses, strains, and mode shape for short-wavelength buckling of 
several different laminates. The compressive stress that corresponds 
to short-wavelength buckling from this analysis for any symmetric 
laminate simplifies to the critical compressive stress for short- 
wavelength buckling from previous studies for unidirectional 
laminates. The equations for the laminate compressive stress 
corresponding to short-wavelength buckling are dominated by matrix 
contributions. The compressive stress corresponding to short- 
wavelength buckling of the outer laminae is half the compressive 
stress that corresponds to short-wavelength buckling of the interior 
laminae. The short-wavelength buckling and delamination of these 
outer laminae characterize the compression failure of [0] g -class 

laminates. The compressive stress that corresponds to short- 
wavelength buckling of some quasi-isotropic laminates is lower than 
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the compressive stress that corresponds to short-wavelength buckling 
of unidirectional laminates with the same number of laminae. This 

O 

lower stress may be due to these quasi-isotropic laminates having 0 
laminae as the outer-most laminae and not having any interior 

O 

0 laminae. The compressive stress that corresponds to short- 
wavelength buckling is a function of fiber volume fraction, and this 
compressive stress is determined from a single curve. All laminates 
in this study buckle into the short-wavelength shear mode. This mode 
shape is dominated by interlaminar shearing, and extensional 
deformations between laminae are negligible. All balanced, symmetric 
laminates have mode shapes with half-waves oriented normal to the 
direction of applied load. Most of the half-wavelengths for short- 

wavelength buckling are on the order of 0(10 1 ). The magnitude of 
these half-wavelengths may influence the experimental compressive 
strengths for composite materials. Some compressive test methods for 
composite materials may inhibit the natural failure mode for a 
laminate by supporting the test specimen along the length in a manner 
that suppresses or prevents short-wavelength buckling. The degree of 
specimen support for different test methods may contribute to the 
scatter in compressive strength data for composite laminates. 

The nonlinear analysis for laminae with short-wavelength 
initial imperfections is used to determine laminate stresses and 
interlaminar strains. This analysis provides the capability to 
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calculate these stresses and strains for a variety of laminates. The 

nonlinear behavior of four laminates with orientations of [0 ] , 

2 s 

[0/90] , [±45], and [+45/0/-45/90] is discussed. The initial 

o o 5 

imperfection for each lamina have the same shape as the laminate's 
short-wavelength buckling mode, and results are presented for 
imperfection-amplitude-to-lamina-thickness ratios of 0.1 and 0.5. The 
compressive-stress verses end-shortening behavior for all laminates is 
similar to the behavior of a wide column. The wide-column response 
for each composite laminate defines a maximum compressive stress that 
corresponds to the maximum load-carrying capacity of the laminate. 

The distribution of the lamina inplane stress resultants, and 

N , across each lamina is discussed. The stress resultant N at 
xy’ x 

the edges of the lamina is greater than at the center of the 

lamina. The stress resultant N is constant across the lamina 

xy 

width. The stress resultant N x is larger than the stress resultant 

N , as expected; however, N is sufficient to initiate failure 

within some laminates. The interlaminar shear strains due to the 
initial imperfections are calculated. The w, displacement 

gradients cause significant interlaminar shear strains ^ KZ ‘ The 

interlaminar shear strains Y are greater than 0.03 for laminate 

compressive loadings that are less than thirty percent of the laminate 
loading for short-wavelength buckling of some laminates. The Y 

Xu 
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for each laminate is largest between the outer two laminae. The 

effects of fiber volume fraction on the nonlinear laminate response 

are presented for [0„] and [±45] laminates. The laminate stiffness 

d s s 

and the maximum Y in the laminate increase as the fiber volume 
xz 

fraction increases. 

A failure criterion for compression-loaded laminates is 

presented. Laminate failures that initiate by outer-lamina buckling, 

by interlaminar shear strains from lamina imperfections, or by inplane 

matrix shearing are included in the criterion. The present linear 

analysis is used to calculate the compressive stress that corresponds 

to outer-lamina buckling, and the present nonlinear analysis is used 

to calculate the compressive stress that corresponds to the 

interlaminar and inplane shear failures. The laminate strength is 

calculated as a function of lamina orientation for [±0] laminates. 

s 

o o 

Compressive failure of [±9] laminates for 0 S 0 < 15 is due to 

s 

outer-lamina buckling when w 0 /t < 0.1 and is due to interlaminar 
shearing when w 0 /t £ 0.1. Compressive failure of these laminates for 
0 0 

15 ^ 9 S 75 is due to inplane matrix shearing. Results from the 

O 

present analysis are not applicable for 0 > 75 . The failure of 

O 

[±0] s laminates with 9 > 75 appears to be initiated by compressive 
failure of the matrix, and this type of failure is not considered in 
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the present model. A simple method called the stiffness-ratio method 

O 

is introduced for predicting the strength of 0 -dominated laminates. 
The stiffness-ratio method uses results from the present nonlinear 
analysis, and the difference between predicted strengths and 
experimental strengths is less than five percent for a limited number 
of laminates. 

4.2 Recommendations for Future Studies 

An analytical study of short-wavelength buckling and shear 
failures in symmetric composite laminates has been described. 
Recommendations for future studies are grouped in three categories: 
analytical extensions; additional experimental verification; and 
applications of the present theory. A useful analytical extension to 
this study would be to generalize the present theory to unsymmetric 
laminates. Such laminates are being considered for use in the next 
generation of commercial transport aircraft structures. Another 
useful analytical extension would be to include material nonlinearity 
in the present theory. Material properties may be a function of the 
compressive load level. These properties are also affected by the 
residual thermal stresses in the laminate and the operating 
temperature of the composite structure. The formulation of the 
present theory is sufficiently general that these analytical 
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extensions could be included by modifying the strain energy 
expressions. 

Additional experimental verification of the present analytical 
results is needed. The effect of the compression test method on a 
laminate's compressive strength and failure mode should be evaluated. 

A study of compression-loaded [±0] laminates might focus on the 

5 

short-wavelength buckling and shear failure mechanisms that were 

0 

described for these laminates. An experimental study of 0 -dominated 
laminates is needed to better evaluate the stiffness-ratio method for 
predicting the compressive strength of these laminates. 

The present theory may be applied in current research. 
Interleaved materials are currently being considered as one of the 
next generation of composite materials. These materials have soft 
adhesive layers that are used (or interleaved) between stiff composite 
layers. Laminates of interleaved materials are similar to the model 
in the present analysis, i.e., alternating soft and stiff layers. The 
present theory may be very useful for predicting the failure 
mechanisms in laminates of interleaved materials. The present theory 
also may be useful for understanding the compressive behavior of woven 
composite laminates. The initial imperfections due to weaving are 
known and may dominate the compressive response of these laminates. 

An investigation of the unique failure mechanisms of woven composite 


laminates needs to be conducted 
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Table 2.1, Typical elastic properties for the graphite-epoxy 
fiber-plate and matrix-foundation. 


Fiber - Plate 

Longitudinal Young's modulus, • E , Msi 

18.50 

Transverse Young's modulus, E ^ , Msi 

1 .64 

Shear modulus, G.^, Msi 

0.87 

Major Poisson's ratio, 

0.30 

Minor Poisson's ratio, v 21 

0.03 f 

Matrix - Foundation [54] 

Young's modulus, E , Msi 

0.59 

Shear modulus, G , Msi 

0.26 

Poisson's ratio, v 

0.36 

ra 



''"calculated using the reciprocal relation 
v 12 = V 21 
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Table 3.1. Short-wavelength buckling results for selected 

two-lamina laminates. 


Laminate 

Orientation 

Compressive 
Stress, ksi 

Compressive 

Strain 

Half- 

Waveleng th, 
in. 

1 — 1 
o 
1 — 1 
to 

289.6 

0.0157 

0.2158 

l>io] s 

283.1 

.0165 

.2188 

l>30] s 

251 .4 

.0342 

.2255 

[+45] s 

253.7 

.0848 

.2066 

[ + 60] s 

279.7 

.1484 

.1544 

1 — 1 
+ 
Co 

0 

1 i 

to 

287.4 

.1741 

.1189 

[90] 

s 

287.5 

.1753 

.1177 









157 


Table 3.2. Short-wavelength buckling results for selected 

four-lamina laminates. 


Laminate 

Orientation 

Compressive 
Stress, ksi 

Compressive 

Strain 

Half- 

Wavelength, 

in. 


^33.5 

0.0234 

0.1905 

co 

i — i 

o 

+1 
1 — » 

433.4 

' 

.0253 

.1910 

[±30] s 

432.7 

.0589 

.1703 

[±45] s 

431 .9 

.1444 

.1450 

[ ± 60] s 

431 .2 

.2289 

.1173 

[±80] s 

430.8 

.2610 

.0969 

C90] s 

430.8 

.2627 

.0973 

[0/90] g 

431 .6 

.0426 

.1336 
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Table 3.3. Short-wavelength buckling results for selected 

eight-lamina laminates. 


Laminate 

Orientation 

Compressive 
Stress, ksi 

Compressive 

Strain 

Half 

Wavelength, 

in. 

co 

r — i 
-=?* 

0 

1 i 

505.4 

0.0273 

0.1854 

[(±45) 2 ] s 

503-7 

.1684 

.1371 

[0/±45/90] 

s 

438.8 

.0589 

.01 40 

[0/+45/90/-45] 

s 

437.2 

.0587 

.0182 

[+45/0/-45/90] 

s 

501 .2 

.0675 

.0477 

[90/+45/0/-45] 

s 

503.6 

.0678 

.1336 
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Table 3.4. Laminate compressive strengths. 


Laminate 

Orientation 

Predicted 
Strength, ksi 

Experimental 
Strength, ksi 

Diff erence^, 
percent 

1 — 1 
o 

L_J 

H 

195.1 

204 

- 4.4 

CO 

i — i 

o 

N 

o 

1 — t 

106.2 

110 

- 3.5 

[0/±45/90] 

s 

64.9 

66 

- 1.7 


(predicted) - (experimental) 
(experimental ) 
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Figure 2.3. Model for a typical N-lamina laminate. 


164 





jj 

« m 

<D J-> 

x c 
n 3 a> 


rH £_ 

as o 


£- t. 


(D 0) 


£» 

CO "H 


e 


i-H CD 
CTJ + 
Q, 

•*H CQ 

o 

C .C 

•rH 4-) 


t- 


a. 5 


■=r 

C\J 

a> 

s- 

3 

b 0 

•r-4 

Du 



165 



Figure 2.5. Short-wavelength buckling mode shapes. 
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Figure 2.6 . Equations (2.69) for the linear analysis of a model with 
anisotropic fiber-plates. 




Figure 3.1. Normalized compressive stress for short-wavelength buckling of 
[0] -class laminate as a function of the number of laminae. 
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Figure 3.2. Normalized laminate compressive stress for short-wavelength 
buckling as a function of fiber volume fraction. 
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Figure 3 « 3 . Mode shape for short-wavelength buckling of a [0] laminate 
(normal waves). 



Figure 3.iJ. Mode shape for short-wavelength buckling of a [ + 10] laminate 
(skewed waves). 
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Figure 3.5. Mode shape for short-wavelength buckling of a [+20] laminate 

o 

(skewed waves) 
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Figure 3 . 6 . Mode shape for short-wavelength buckling of a [+30] a laminate 
(skewed waves). 
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Figure 3.7. Mode shape for short-wavelength buckling of a [+45] s laminate 
(skewed waves). 
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Figure 3.9. Mode shape for short-wavelength buckling of a [+80] s laminate 
(skewed waves). 
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Figure 3.10. Skew angle for the laminate mode shape as a function of fiber 
orientation. 
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Figure 3.11. Compressive strain for short-wavelength buckling as a function of 
laminate half-wavelength. 
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Figure 3.12. Compressive strain for short-wavelength buckling as a function of 
laminate half-wavelength for a [0 ] laminate (expanded ordinate). 




Figure 3.13. Half-wavelength of the buckling mode as a function of lamina 
orientation for a [+e] laminate. 



Figure 3.1^. Half-wavelength of the buckling mode as a function of lamina 
orientation for a [±e] laminate. 




Figure 3.15. Half-wavelength for short-wavelength buckling as a function of 
laminate fiber volume fraction. 
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Figure 3.16. Normalized compressive stress versus end shortening for 
laminate. 
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Figure 3.17. Stress resultant for a 0 lamina in a [0 ] laminate (w 0 /t 
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Figure 3.18. Distribution of the w displacement for the outer-most 0 lamina in 
a CO] laminate along x = A/2 (w 0 /t = 0.5). 
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Normalized maximum w displacement, 
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Figure 3-20. Normalized compressive stress versus end shortening for a [0/90] 
laminate. 
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Figure 3.22. Normalized compressive stress versus end shortening for a [±4&] 
laminate. 
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Figure 3.23. Stress resultants for a H5 lamina in a [±45] , laminate 
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Figure 3.2^. Distribution of the w displacement for the outer-most 45 lamina 
in a [±45] laminate along x = A/2 (w o /t = 0.5). 



192 



Figure 3.25. Maximum w displacement as a function of end shortening for 
[±^5] laminate. 
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[+45/0/-45/90] laminate. 
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laminate for u 0 /(u 0 ) i = 0.2b (w 0 /t 
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Figure 3.28. Stress resultants for 0 and 45 laminae in a [+45/0/-45/90] 
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Figure 3.31. Interlaminar shear strains in a [0 o ] laminate for 
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Figure 3.3^. Interlaminar shear strain at the 0 /90 interface along x = 0. for 

[0/90] laminates for u 0 /(u D ) = 0.20. 

s cr* 
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Figure 3.35. Interlaminar shear strain_at t.he +45 /-45 interface along 
for [±45] g laminates for u o /(u o ) cr = 0.20. 
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Figure 3.37. Effects of fiber volume fraction on the normalized compressive- 
stress - end-shortening behavior for a [0 ] laminate. 
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Figure 3.40. Effects of fiber volume fraction on the normalized compressive- 
stress - end-shortening behavior for a [+45] laminate. 
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a function of the end shortening for a [±i)5] laminate. 
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Figure 3. **3. Compressive strength for a [±0], laminate as a function of lamina 
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